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1.  Background 


This  is  the  final  report  of  the  project  "Study  of 
Mark's  Electrothermodynamic  (ETD)  Power  Generator",  spon¬ 
sored  by  the  Department  of  Energy  (DOE) .  The  primary  ob¬ 
jective  of  this  effort  has  been  to  provide  a  detailed 
theoretical  analysis  of  a  type  of  electrohydrodynamic  (EHD) 
power  generator  proposed  by  the  Marks  Polarized  Corporation. 
Once  a  sound  analytical  basis  was  developed,  it  was  possible 
to  ascertain  the  limitations  of  the  proposed  system  and  to 
seek  improvements. 

The  results  of  our  first  study  are  reported  in  detail 
in  Ref.  (1)  and  summarized  more  concisely  in  Ref.  (2) . 

That  work  indicated  the  originally  proposed  design  would 
yield  unacceptably  poor  performance  and  suggested  a  certain 
design  revision.  Further  work  was  done  to  analyze  such  a 
design  revision  and  is  reported  in  the  Appendix  to  Ref.  (1) . 
An  improved  simplified  and  shortened  version  of  that  work  is 
reported  in  Ref.  (3).  Unfortunately,  the  new  results  indi¬ 
cated  that  the  revised  design,  while  showing  somewhat  im¬ 
proved  performance,  still  would  not  appear  sufficiently 
competitive  in  overall  efficiency. 

In  the  first  study,  some  minor  effects  of  compressibili¬ 
ty  of  the  medium  were  neglected.  This  was  done  in  order  to 
simplify  the  calculations.  Also,  it  has  been  conjectured  by 
Marks  and  others  that  it  might  be  possible  to  improve  greatly 
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the  performance  of  EHD-generators  which  utilizes  an  ejector 
by  employing  a  two-fluid  cycle.  Finally,  performance  gains 
that  a  diffuser  might  bring  about  appeared  worthy  of  study. 

To  this  end,  a  second  study  for  a  two-fluid  cycle  was  under¬ 
taken  in  which  a  diffuser  is  introduced  just  prior  to  the 
EHD  generator  and  where  full  account  is  taken  of  the  com¬ 
pressibility  of  both  fluids  at  all  locations. 

The  principal  earlier  effort  along  the  lines  of  our 
second  study  is  Ref.  (4),  by  Huberman,  et.  al ,  but  our 
present  work  has  gone  beyond  anything  undertaken  in  Ref. 

(4).  Although  a  two-fluid  generator  is  analyzed  in  Ref. 

(4) ,  including  the  compressibility  of  both  fluids,  Huberman 
et.  al,  make  no  attempt  to  study  the  complete  thermodynamic 
cycle.  There  are,  moreover,  a  number  of  questionable 
assumptions  which  cast  serious  doubt  in  the  validity  of  the 
results  given  in  Ref.  (4)  . 

The  computational  scheme  that  we  have  developed  was 
reported  in  Ref.  (5).  There  a  methodical  derivation  of 
the  working  formulas  is  given  together  with  some  preliminary 
results.  Some  further  results  are  given  in  Ref.  (6)  . 

The  present  report  updates  the  material  given  in  Refs 

(5)  and  (6) .  In  particular,  we  have  found  that  instead  of 
writing  an  energy  balance  across  the  condenser,  it  is  more 
convenient  to  use  a  momentum  balance  as  shown  in  Section  II 
of  this  report.  We  have  reproduced  herein  the  complete  analy¬ 
sis  in  its  latest  version. 
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A  sample  calculation  for  the  two-fluid  system  is  also 
included  together  with  the  computer  program  listing. 

Figure  1.1  is  a  schematic  of  the  generator  configura¬ 
tion  studied.  In  Marks'  design,  the  centerline  of  the 
channel,  which  is  shown  in  Fig.  1.1  as  a  straight  line,  is 
curved  sc  as  to  form  a  closed  loop. 


2. 


Introduction 


For  purposes  of  analysis,  the  overall  thermodynamic  pro¬ 
cess  in  the  two-fluid  EHD  generator  may  be  represented  by  the 
schematic  flow  diagram  shown  in  Fig.  2.1. 

Primary  fluid  enters  the  ejector  at  station  1,  secondary 
fluid  enters  at  station  2  and  the  mixture  leaves  at  station  3. 

The  mixture  then  passes  through  a  diffuser,  process  3-*-4 , 
the  basic  purpose  of  which  is  to  reduce  the  kinetic  energy  of 
the  flow  and  thereby  decrease  the  associated  friction  losses 
downstream  of  the  diffuser.  On  the  other  hand  the  diffuser 
introduces  certain  losses  of  its  own  so  that  some  care  is 
necessary  to  optimize  the  design  in  this  respect. 

The  gas  mixture  then  passes  through  the  electrical  power 
section,  process  4-*-5,  from  which  the  gross  electrical  power  is 
extracted. 

Finally,  the  gas  mixture  enters  the  condenser/separator 
at  station  5.  The  primary  fluid  is  separated  from  the  gas 
mixture  by  condensation  and  leaves  at  station  6.  The  fluid  at 
this  point  is  assumed  to  be  compressed  liquid  at  a  known  static 
pressure  Pg  and  a  known  static  temperature  Tg  .  Kinetic 
energy  at  station  6  is  regarded  as  negligible.  The  secondary 
fluid  leaves  the  condenser/separator  and  enters  the  ejector  at 
station  2.  We  assume,  moreover,  that 

P6  =  P2  =  Px  (2.1) 

and 

T6  =  T2  (2.2) 
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|.2.l  SCHEMATIC  FLOW  DIAGRAM  OF  TWO  -  FLUID  EHD  CYCLE. 


The  primary  fluid  which  leaves  the  condenser/separator 
at  station  6  is  circulated  back  through  the  pump  and  through 
the  boiler/superheater,  process  7-+-1. 

It  should  also  be  explained  that  the  primary  fluid  at 
station  1  is  a  condensible  vapor  of  high  molecular  weight 
which  contains  fine  liquid  droplets.  Unless  specifically 
stated  otherwise  the  quality  z  at  this  location  is  taken 
as  0.95.  The  important  thermodynamic  static  properties  at 


station  1  can  be  expressed  in  the  following  form. 

h^  =  hf  +  z  hfg  Enthalpy 

S1  =  sf  +  2  sfg  Entropy 


(2.3) 


(2.4) 


Here  subscript  f  denotes  the  saturated  liquid  and  subscript 
fg  denotes  the  change  between  saturated  liquid  and  saturated 
vapor.  All  properties  are  evaluated  at  a  specified  pressure 
P.^  ;  the  corresponding  saturation  temperature  T^  is  then 
also  known. 

The  corresponding  stagnation  properties  Tg^  and  hs^ 
are  fixed  when  the  stagnation  pressure  P  ^  is  specified. 

Of  course  the  entropy  has  the  same  value  s^  at  the  stagna¬ 
tion  condition  as  at  the  corresponding  static  state.  Thus 
specification  of  Pgl  and  s^  suffices  to  fix  all  other 
stagnation  properties. 

While  it  is  convenient  to  start  the  analyses  by  stipula¬ 
ting  the  quality  at  station  1,  it  is  also  advisable  later  in 
the  calculation  to  determine  the  corresponding  quality  at 


station  4,  the  entrance  to  the  electrical  working  section. 

It  is  essential  that  the  fluid  entering  the  working  section 
contain  finely  dispersed  liquid  droplets  of  the  proper  size. 
These  droplets  carry  the  electrical  charges  that  are  the 
central  feature  of  EHD  power  generation.  We  assume  tentatively 
that  acceptable  droplet  characteristics  are  obtained  by 
maintaining  the  quality  at  the  inlet  to  the  working  section 
in  the  range  0.92  to  0.98.  These  limits  are  only  estimates. 

The  thermodynamic  analysis  would  be  greatly  simplified 
if  both  fluids  were  perfect  gases  with  constant  specific 
heats.  In  fact  only  the  secondary  fluid  satisfies  this  re¬ 
quirement.  Nevertheless,  we  can  define  an  "equivalent  perfect 
gas"  which  adequately  approximates  the  essential  thermodynamic 
properties  of  the  primary  vapor/liquid  mixture.  When  carried 
out  judiciously,  this  procedure  simplifies  the  analysis  at  the 
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cost  of  only  a  small  loss  in  numerical  accuracy  . 

The  ratio  of  specific  heats  y  of  the  "equivalent  perfect 
gas"  may  be  deduced  from  the  important  relation 


<Psl/Pl>  *  <Tsl'Tl> 


Y1/(Y1  -  1) 


(2.5) 


Solving  this  for  y^  gives 


to(Tsl/Tl) 

Yi  '  Wp^Tp^I 


(2.6) 


Another  fundamental  property  of  the  "equivalent  perfect 
gas"  is  its  specific  heat  Cpl  which  may  be  evaluated  from 
the  expression 


(2.7) 
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'Pi 


{hsl  “  hl> 
(Tsl  -  Tl} 


where  all  quantities  on  the  right  are  now  known. 

It  then  follows  from  standard  perfect  gas  relations  that 
the  gas  constant  and  the  molecular  weight  W1  of  the 

"equivalent  perfect  gas"  are,  respectively. 


R1  " 


'Yi  -  1\ 


Si 


(2.8) 


and 


w  =  5 

W1  R, 


(2.9) 


where  R  =  universal  gas  constant 
=  8317  joules/kg-mole  °K 

Eqs .  (2.6)  through  (2.9)  complete  the  definition  of  the 
"equivalent  perfect  gas"  which  adequately  simulates  the  essen¬ 
tial  thermodynamic  properties  of  the  primary  vapor/liquid 
mixture. 

It  should  also  be  explained  that  the  secondary  fluid  at 
station  2  is  predominantly  a  noncondensing  gas  of  low  molecu¬ 
lar  weight  primary  vapor. 

It  can  be  shown  that  the  mass  ratio  oi  of  condensible 
primary  vapor  to  dry  secondary  gas  at  station  2  is  given  by 
the  expression 


01 


WG  TP 


2 


vv 

-  W  J 


(2.10) 
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where 


mu  iwy 


Wv  =  true  molecular  weight  of  condensible  vapor  (not  to 
be  confused  with  the  "equivalent  molecular  weight 
W^"  considered  earlier) 

WG  =  molecular  weight  of  noncondensing  gas 
Pv(Tg)  =  vapor  pressure  of  condensible  primary  fluid  at  known 
temperature  Tg 

P2  =  known  static  pressure  at  station  2 
The  present  analysis  is  restricted  to  conditions  under 
which  <u  is  very  small  compared  with  unity  and  may  be  neglec¬ 
ted.  Thus  the  fluid  at  station  2  may  be  treated  as  dry  gas. 
This  greatly  simplifies  the  analysis.  Eq.  (2.10)  is  useful  for 
verifying  that  this  assumption  is  indeed  satisfied  in  any  par¬ 
ticular  instance. 

It  should  be  added  that  if  it  later  becomes  of  interest 
to  make  calculations  for  circumstances  under  which  parameter 
on  is  not  negligible,  the  present  analysis  can  in  fact  be 
generalized  to  include  this  effect.  For  the  present,  however, 
we  prefer  to  deal  only  with  the  simpler  situation  in  which  co 
may  be  neglected. 
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3.  Ejector 

The  ejector  receives  a  primary  stream  of  high  molecular 
weight  gas  at  station  1,  a  secondary  stream  of  low  molecular 
weight  gas  at  station  2  and  discharges  the  resulting  mixture 
at  station  3.  Static  conditions  at  these  three  stations  are 
designated  by  subscripts  1  ,  2  ,  3  .  Stagnation  conditions 

at  the  corresponding  stations  are  designated  by  subscripts 
si  ,  s2  and  s3  . 

The  following  quantities  are  arbitrarily  specified  or 
known  at  stations  1  and/or  2,  namely, 

Molecular  weights :  and  5*2 

Ratios  of  specific  heats:  y^  and  y2 
Stagnation  pressure  of  primary  jet:  Pgl 
Stagnation  temperature  of  primary  jet:  Tgl 
Static  pressures:  P1  *  P2  (=  Pg) 

Static  temperature  of  primary  jet:  T^  (=  saturation  tempera¬ 
ture  of  primary  fluid  at  pressure  P^) 

Static  temperature  of  secondary  jet:  T2  (  =  Tg) 

Mass  flow  ratio:  =  m^/Cif^  +  ra2)  =  x 

Velocity  ratio:  V^/V2  =  ¥ 

We  treat  the  fluids  at  stations  1,  2,  3,  4  and  5  as  per¬ 
fect  gases  with  constant  specific  heats. 

The  velocity  ratio  of  the  ejector  may  be  developed  as 
follows: 


al 

Ml  * 

ft 

R1 

Tsl 

M;  -  _i 

a2 

M2  1 

'  Y2 

r2 

T2  1 

J*77h  ’  M2 

(3.1) 
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We  can  now  use  this  result  and  other  standard  perfect 
gas  relations  to  complete  the  calculation  of  key  thermody¬ 
namic  properties  at  stations  1  and  2.  It  is  convenient  to 
arrange  these  calculations  sequentially  as  follows : 


M 


1 


(3.2) 


M 


2 


fV2\J~Yl  R1  *sl  M1  -  - 

Kj'l  Y2  R2  T2 


T2  (ts2/t2) 
^2 


(3.3) 

(3.4) 

(3.5) 

(3.6) 

(  3.7) 


Before  we  can  determine  the  actual  conditions  at  station 
3,  the  outlet  of  the  ejector,  it  is  first  necessary  to  analyze 
two  other  hypothetical  cases  as  indicated  schematically  in 
Fig.  3.1.  In  all  three  of  these  cases  the  flow  is  taken  as 
steady  and  adiabatic. 

The  first  case,  that  shown  in  Fig.  3.1(a) ,  involves  a 
hypothetical  device  which  receives  two  gas  streams  at  stations 
1  and  2  and  discharges  the  resulting  mixture  at  station  x  .  In¬ 
let  conditions  at  1  and  2  are  identical  to  those  of  the  actual 
ejector. 
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Hypothetical  Minimum 
Loss  Device 
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Fig.  3.1  COMPARISON  OF  ACTUAL  EJECTOR  WITH  IDEAL 
EJECTOR  ANO  WITH  HYPOTHETICAL  MINIMUM 
LOSS  DEVICE. 


Recall  that  the  adiabatic  mixing  of  two  different  gases  is 
an  inherently  irreversible  process  that  always  involves  a  cor¬ 
responding  entropy  increase.  We  stipulate  that  the  first  case 
is  such  that  the  only  irreversibility  which  occurs  is  that  asso¬ 
ciated  with  this  mixing.  We  term  this  mass  mixing  to  distin¬ 
guish  it  from  another  type  of  mixing  considered  below  which  we 
term  momentum  mixing. 

The  second  case,  that  shown  in  Fig.  3.1(b),  involves  an 
ideal  ejector  which  receives  two  gas  streams  at  stations  1  and 
2  and  discharges  the  resulting  mixture  at  station  y.  Again, 
the  inlet  conditions  at  1  and  2  are  identical  to  those  of  the 
actual  ejector.  This  case  satisfies,  among  other  relations,  the 
idealised  one -dimensional  momentum  equation  for  frictionless, 
constant  area  flow.  It  can  be  shown  that  in  this  case  there  is 
an  overall  entropy  increase  which  includes  not  only  the  previ¬ 
ously  mentioned  effect  of  mass  mixing  but  also  a  further  increase 
associated  with  momentum  mixing. 

The  third  case,  that  shown  in  Fig.  3..1(c),  represents  the 
actual  ejector  itself  which  receives  input  streams  at  stations 
1  and  2  and  which  discharges  the  resulting  mixture  at  station  3. 

All  three  of  the  above  cases  refer  to  constant  area  pro¬ 
cesses  in  the  sense  that 

Ax  "  Ay  =  A3  ■  (Ax  +  A,)  (  3.8) 

Once  the  hypothetical  operating  conditions  at  stations  x 
and  y  have  been  found,  it  then  becomes  possible  to  express  the 
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corresponding  actual  conditions  at  station  3.  This  is  accom¬ 
plished  through  the  use  of  an  ejector  effectiveness  nE  . 

This  parameter  is  defined  more  explicitly  in  the  later  analysis. 
Values  of  n£  must  be  estimated  by  reference  to  test  data  on 
ejectors  (i.e.,  Ref.  8/  9 •  10). 
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4.  Adiabatic  Mixing  of  Two  Streams 

The  process  in  the  ejector  is  treated  as  the  adiabatic 
mixing  of  two  streams.  The  following  relations  can  be  shown 
to  apply. 


R  =  8315  Joule/kg -mole  °K 


r_  =  R_ 

R 1  wx 


R  -  R 
2  "  ^ 


R„  =  R  =  R,  =  X  R.  +  (1  -  x)  R. 
X  y  3  1  * 


Y1  R1 


'pi  (Y-,  -  l) 


y2  R2 

Cp2  "  (Y,  -  1) 


:  =  c  =c,  =  xc1  +  (i-x)c_ 

px  py  p3  pi  p2 


IE2. 


Y*  *  Yy  Y3  -  lcp3  -  R3T 


P  -  T  =  T  -  =  ■==-  [x  C  .  T  .  +  (1 
sx  sy  s3  Cp3  pi  si 


Ax  =  Ay  =  A3  =  (A1  +  A2} 


(4.1) 

(4.2) 

(  4.3) 


(4-4) 


(  4.5) 


(  4.6) 

(4.7) 


x) 


Ts2  ^ 


(4.8) 

(4.9) 

(4.10) 


Notice  that  the  quantities  R  /  C  ,  y  »  T  and  A  all 

p  s 

have  the  same  values  at  stations  x,  y  and  3. 

This  fact  may  be  used  to  simplify  many  of  the  subsequent 
expressions.  It  is  always  permissible  to  substitute  subscript 
3  for  subscripts  x  or  y  on  any  of  these  quantities. 
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5.  Mass  Flow 


It  is  convenient  temporarily  to  omit  station  subscripts 
1,  2,  3,  x,  y  and  thereby  develop  certain  needed  mass  flow 
relations  in  generalized  form.  Thus  the  mass  flow  across  an 
arbitrary  station  may  be  written  and  developed  as  follows. 

m  =  pav  =  A  (/Trt  m)  =  PAM  JjL 


P  .  A  T  /T 
s  ^ 

/ rts7y  "(ps7p) 


This  may  be  rewritten  as 


P  A 


FV7 


where  the  auxiliary  function  f(M)  is  defined  as  follows. 


f(M)  = 


>/Ts/T 

TP^PT 


M  =  M 


(Y  +  1) 
(y  -  1) 


Eq.  (5.2)  can  be  applied  specifically  to  stations  1,  2, 
x,  y,  3  as  follows 

P  .  A. 

m,  =  s-  1  f  (M. )  (5. 

/r  t  77"  x 

/K1  sl/Yl 


Ps2  ^2 

*5  *  f  (M_) 


'2  s2'  '2 


All  quantities  on  the  right  side  of  Eq.  (5.13)  except 
P  are  now  known.  The  method  of  determining  P  is  ex- 

S  X  sx 

plained  in  a  later  section;  refer  to  Eq.  (5.15) .  Once  P 

S  X 

has  been  found,  Eq.  (5.13)  fixes  f(M^)  whereupon  Eq.  (5.3) 
fixes  M  itself.  The  solution  of  Eq.  (5.3)  for  M  when 
f(M)  is  known  involves  an  iterative  procedure  which  is  ex¬ 
plained  later  in  this  section. 

From  Eqs.  (5.6),  (5.7),  (5.8)  and  (5.10)  we  also  find  that 

P  f (M  )  =  P  f (M  )  =  P  ,  f (M,)  (5.14) 

sx  x  sy  y  s3  3 


Assuming  P  and  M  known,  this  relation  along  with  Eq.  (5.3) 

SX  X 

fixes  P  when  M  is  specified;  it  also  fixes  f(M_,)  and 
sy  y  3 

when  Pg2  is  specified.  The  methods  of  calculating 
and  Pg3  are  explained  later;  refer  to  Eqs.  (7.8)  and  (8.14). 

Fig.  5 . 1  is  a  rough  sketch  of  the  function  f(M)  as  defined 
by  Eq.  (5.3).  This  function  vanishes  at  M  =  0  and  at  M  =  »  . 
The  maximum  value  occurs  at  M  =  1  and  equals 


(5  .15) 


It  is  clear  from  the  figure  that  in  the  range  0<f<  fmax  • 
Eq.  (5.3)  has  two  positive  real  roots,  one  subsonic  and  the 
other  supersonic.  Usually  the  supersonic  root  must  be  discarded 
because  it  is  not  consistent  with  a  further  constraint  imposed 
by  the  second  law  of  thermodynamics. 
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In  order  to  solve  Eq.  (5.3)  for  M  when  values  of  f 
and  y  are  specified,  we  utilize  an  iterative  procedure  based 


on  Newton's  method.  Let  M  denote  the  nth  trial  value  of 

n 

M  .  We  then  define  a  function  F(M  )  and  its  derivative 

n 

I 

F  (M  )  as  follows, 
n 


F(Mn)  =  f 


1  +  Y -  M  “ 
2  n 


(Y  +  1) 

2Ty  -TV 


-  M 


n 


(5.16) 


(H-i)  «» 

f1  +  (Y  2  *)  Mn2  J 

(3  -  Y) 


-  1  (5.17) 


The  (n  +  l)*t  trial  value  of  M  may  now  be  taken  as 

F  (M  ) 


M 


(n  +  1)  ~  Mn  ‘  ^ 


n 


F  (Mn> 


(5.18) 


The  cycle  of  calculations  defined  by  Eqs.  (5.16)  ,  (5.17) 
and  (5.18)  is  repeated  until  the  result  for  converges  to 
a  stable  value  at  the  desired  level  of  accuracy. 

Thus 


M  *  Aim  M 
n 

n  -►  00 


(5.19) 


It  is  readily  apparent  that  the  value  of  M  so  obtained 
satisfies  Eq.  (5.3)  for  the  prescribed  values  of  f  and  y  . 
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6 .  Entropy 

Consider  the  case  shown  in  Fig.  3.1(a).  The  specific 
entropy  sx  of  the  discharged  gas  may  be  expressed  in  the 
form 


where  Px^  and  Px2  are  the  partial  pressures  of  the  two  com¬ 
ponents  of  the  gas  mixture.  Symbols  PQ  and  Tq  denote  the 
pressure  and  temperature  of  the  ambient  atmosphere.  The  en¬ 
tropy  of  each  component  is  assigned  the  value  zero  at  this 

reference  state  P„  ,  T 

o  o 

Incidentally,  the  primary  fluid,  being  condensible,  might 
well  exist  only  in  the  liquid  state  at  the  reference  condition 
PQ  ,  Tq  .  Nevertheless,  so  long  as  it  is  in  the  gaseous  state 
at  station  x  ,  Eq.  (6.1)  may  still  be  used. 

The  partial  pressures  of  the  two  components  are  propor¬ 
tional  to  the  respective  mol  fractions.  Hence 


(6  .2) 


(6  .3) 
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4 


i  — 


Upon  substituting  Eqs.  (6.2)  and  (6.3)  into  (6.1)  and 


simplifying  we  may  obtain  the  result  in  the  form 


SX  "  3nt  +  Cp3  “rh  Mp- 


(6.4) 


where 

In  R3  -  x  R1  Hn  (x  R^)  -  (  1  -  x)R2  &n  [(1  -  x)R2] 

(6.5) 


It  can  be  shown  that  this  quantity  sm  represents  the 
entropy  increase  associated  with  the  mass  mixing  of  the  two 
different  gases. 

Moreover,  since  the  entropy  of  the  stagnation  state  s^ 
is  by  definition  identical  to  that  of  the  corresponding  static 
state  x  ,  we  may  replace  T  and  in  Eq.  (6.4)  by  T 

X  X  3* 

and  P  ,  respectively.  Thus  we  obtain 

SX 

/p  \ 

3  <6-6> 


In  previous  studies  of  this  series,  we  have  considered 
only  the  special  case  in  which  the  same  fluid  is  employed  for 
both  the  primary  and  secondary  streams.  In  that  case  we  have 

Rx  »  R2  =  R3  =  R  (6.7) 


Cpl  “  Cp2  ~  Cp3  =  Cp 


(  6.8) 


Y 


1 


=  Y2  =  Y3  *  Y 


(  6.9) 


Moreover,  in  this  special  case  there  is  no  mass  mixing,  so 
that  we  must  set 
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0 


(6-10) 


s 

m 


Consequently,  Eq.  (6.6)  now  simplifies  to 


s  =  0  + 
x 


(6-11) 


Next  reverting  to  the  case  of  two  different  gases,  let 
At  denote  a  time  interval  during  which  unit  mass  crosses 
station  x  .  During  this  same  interval  x  units  of  mass  cross 
station  1  and  (1  -  x)  units  of  mass  cross  station  2.  The 
entropy  which  leaves  across  station  x  must  equal  the  entropy 
which  enters  across  stations  1  and  2  plus  the  entropy  increase 
caused  by  the  mass  mixing.  Thus 


3x  *  +  *  S1  +  (1  -  x)s2 


t  6.12) 


where 


S1  =  Cpl  ln 


s2  ~  Cp2 


_  p 

in  | 

Tc  ) 

Ri 

lPo  . 

—  P 

in 

(— 

yr) 

R2 

1-0 

(6.13) 


(  6.14) 


Upon  eliminating  between  Eqs.  (  6.6)  and  ('6.12)  ,  we 

find  that  s_  also  cancels  from  the  result.  It  is  then  a 
m 

simple  matter  to  solve  for  the  exit  stagnation  pressure  in  the 
form 


(  6.15) 
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m 


Thus  the  sequence  defined  by  Eqs.  (6.13),  (6.14)  and 

(6.15)  now  fixes  P  .  Next  referring  back  to  Eq.  (5.13), 

S  X 

we  can  calculate  the  value  of  f  (M^)  ;  finally,  from  Eqs. 

(5.16) ,  (5.17),  (  518)  we  can  calculate  itself.  Once 

P  ,  T  =  T  ,  and  M  are  known,  it  is  a  simple  matter  to 

sx  '  sx  s3  x 

calculate  the  corresponding  conditions  Px  ,  from  the 

usual  perfect  gas  relations. 
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p 


g(M)  =  - % — «-  A  +  (X_lJ:)M2  (7.5) 

(1  +  YM^)  Z 

With  this  notation  Eq.  (7.3)  may  be  rewritten  in  the 

form 


x  /R1  Tsl  (1  -  x)  /*2  Ts2  _  1  7*3  Ts3 

gTHp"/  yl  g(M2)  /  y2  g(My)  /  y3 


(  7.6) 


Rearramging  gives 


g<M  >  =  / 


g(Mi} 


T 

1  Si 


(1  -  X) 

g(M2V 


-l 


(  7  7) 


Since  all  quantities  on  the  right  are  known,  Eq.  (7.7) 
fixes  g(My)  .  Then  itself  can  be  found  by  inverting 

Eq.  (7.5).  Fortunately,  an  explicit  solution  is  possible  in 
this  case  as  Eq.  (7.5)  can  be  reduced  to  a  quadratic  in  M2  . 
The  result  is 


M2  =  (1  -  2y?2)  +  A  V  *  UZ  (7.8) 

1  -  Y ( 1  -  2Yg  ) 

The  general  character  of  the  function  g(M)  is  sketched 

in  Fig.  7.1.  Notice  that  g(M)  vanishes  at  M  =  0  ,  reaches 

its  peak  value  9mx  at  M  =  1  ,  and  decreases  toward  the 

limit  g  as  M**»  .  It  cam  be  shown  that 
^00 


gmax 


1 

✓2(y  +  1) 


(7.9) 


g  = 


(7.10) 


27 


Fig.  71  APPROXIMATE  SKETCH  OF  FUNCTION  g(M) 


Fig.  7.1  discloses  the  existence  of  a  supersonic  root  over 
the  range  <?oo<<?<<3max’  This  root  should  not  be  accepted  until 
a  check  calculation  is  made  to  determine  whether  the  result  is 
consistent  with  the  second  law  of  thermodynamics.  The  subsonic 
root  is  found  by  retaining  only  the  negative  sign  before  the 
radical  in  Eq.  (7.8). 

Once  M  has  been  found  from  Eq.  (7.8),  the  corresponding 
value  of  P  may  be  found  from  Eqs.  (5.3)  and  (5.14).  This 
suffices  to  fix  all  properties  at  states  y  and  sy  . 


29 


8.  Availability 

The  entropies  and  S£  of  the  two  input  streams 

were  defined  earlier  in  Eqs.  (6.13)  and  (6.14).  The  corres¬ 
ponding  steady  flow  availability  functions  with  respect  to  an 


ambient  atmosphere 

be  written 

at  pressure 

P 

o 

and  temperature 

Tq  may 

h 

=  C  . 
pl 

(Tsl  - 

V 

-  T  s, 
o  1 

(8.1) 

*2 

=  C  -5 

P2 

(Ts2  * 

V 

-  T  S, 
o  1 

(8.2) 

Consequently  the  total  available  energy  entering  the  sys¬ 
tem  becomes 


■J>m  "  x  ^  +  (1  “  x)^2 


(  8.3) 


The  respective  entropies  of  the  streams  leaving  at  stations 
x  ,  y  and  3  are 


s  *  s  +  C  ,  in  -  R3  Jin  (  *25 

*  m  p3  \To  /  3  \Po  y 


(  8.4) 


sy  ’ s*  +  S3 r2  -  e3  *» 


(8.5) 


S3  *  +  =P3  «»  fe3)  -  *3 


(8-6) 


where  sm  has  been  previously  defined  in  Eq.  (6.5). 

The  corresponding  availabilities  of  the  streams  leaving 
at  stations  x,  y  and  3  may  therefore  be  written 
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(8.7) 


p  =  c  , 

*x  p3 

{Ts3  - 

V 

-  T 

O 

s 

X 

(8.7) 

*’y  =  cp3 

(Ts3  - 

V 

-  To 

s 

y 

(8.8) 

*3  *  Cp3 

(Ts3  - 

V 

-  T 

o 

S3 

(8.9) 

Now  consider  the  losses  of  availability  listed  below. 

By  substituting  the  above  expressions  for  the  quantities  on 
the  left  and  simplifying  we  readily  obtain  the  expressions 
shown  on  the  right.  Thus 

<*m  -  V  *  To  sm  < 

<*x  '  V  *  R3  To  l8-U) 

<*x  -  *3>  "  a3  To  “ 

These  results  are  very  significant.  Eq.  (8.10)  defines 
the  loss  of  available  energy  caused  by  the  mass  mixing  of  the 
two  different  gases.  Eq.  (8.11)  defines  the  further  loss  of 
available  energy  caused  by  the  momentum  mixing  in  an  idealized 
ejector.  Eq.  (8.12)  defines  the  augmented  loss  of  available 
energy  caused  by  the  momentum  mixing  in  the  real  ejector. 

It  is  useful  to  postulate  an  empirical  relation  between 
these  last  two  losses  which  is  shown  below  on  the  left.  The 
expression  on  the  right  then  follows  from  Eqs.  (8.11)  and 
(8.12).  Thus 


sx 


(  8.12) 


(’p  -  p  ) 

T  V  XT 


In  (P  /P  ) 

otr  a  v 


We  term  nE  the  ejector  effectiveness.  Its  value, 
which  is  always  less  than  unity,  must  be  estimated  from  appro¬ 
priate  test  data  on  ejectors. 

If  we  treat  nE  as  known,  Eq.  (8.13)  fixes  .  The 

solution  is  simply 


We  cam  next  find  f(M^)  from  Eq.  (5.14)  and  from 

Eqs.  (5.16),  (5.17)  and  (5.18).  This  suffices  to  fix  all  pro¬ 
perties  at  states  3  and  s3  thereby  completing  the  analysis  of 
the  ejector. 
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Diffuser 


9  . 

It  is  convenient  to  analyze  the  diffuser  in  terms  of  the 

concept  of  availability  in  steady  flow.  However,  we  wish  to 

deal  with  the  characteristics  of  the  diffuser  itself  and  these 

cannot  depend  on  the  arbitrary  values  PQ  ,  Tq  which  happen 

to  characterize  the  condition  of  the  ambient  atmosphere.  We 

can  achieve  our  objective  by  defining  the  availability  with 

respect  to  a  reference  state  ,  T3  which  characterizes 

the  diffuser  rather  than  with  respect  to  the  ambient  state 

P  ,  T  . 
o  '  o 

Under  these  circumstances  we  may  write  the  appropriate 
availabilities  at  diffuser  inlet  and  outlet  as  follows 


'^s3  “  Cp3 


(Ts3  -  t3) 


(9.1) 


s4 


=  C 


P3 


(T 


s3 


-  T3} 


(9.2) 


These  equations  make  use  of  the  fact  that  Tg4  =  T  ^  • 
Also,  states  s3  and  3  are  at  the  same  entropy  so  that 


(ss3  -  s3> 


(9.3) 


From  Eqs.  (9.1)  and  (9.3)  we  conclude  that 


This  result  shows  that  the  available  energy  at  the 


diffuser  inlet  is  simply  the  inlet  kinetic  energy  itself. 

The  loss  of  availability  through  the  diffuser  may  be 
found  by  subtracting  Eq.  (  Sl2)  from  (9.1).  The  result  is 


(*s4  "  W  =  W  =  +  R3  T3  in  {  F 


s  3 
s4 


We  now  define  the  diffuser  effectiveness  as 


s4 

PS3. 


*  n, 


so  that 


»S3  -  ***\ .  a_ . 
*33  /  *s3 

This  reduces  readily  to 


(l  -  nD> 


r3  t 


Solving  for  the  pressure  ratio  gives 


e 


(1  -  nD) 


(9.5) 


(  9.6) 


(  9.7) 


(  9.8) 


(  9.9) 


This  is  the  result  required.  It  fixes  Pg4  when  Pg3  , 
and  nD  are  specified. 
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Assuming  the  area  ratio  A3/A4  to  be  arbitrarily  spec¬ 
ified,  we  can  find  from  the  continuity  relation.  Refer¬ 

ring  to  Eq.  (5.8)  we  may  write 


m3  * 


Ps3  A3 

/*3  ^s7*3 


f(M3)  = 


PS4  A4 
^3  Ts3^3 


f(M4)  =  A4 


(9.10) 


Consequently 


f(M4) 


(9.11) 


Eq.  (9.11)  fixes  f(M4)  .  Then  M4  follows  in  the  usual 
way  from  Eqs.  (5.16),  (5.17)  and  (5.18).  The  result  suffices 
to  fix  all  properties  at  states  4  and  s4  .  Thus 


1  + 


"  Ps4/(Ps4/P4) 


(9-12) 

(9.13) 

(9.14) 

(9.15) 


The  configurations  studied  in  earlier  work  did  not  incor¬ 
porate  a  diffuser.  For  the  purpose  of  comparing  the  results 
of  the  present  analysis  with  corresponding  earlier  work,  it  is 
desirable  to  be  able  to  eliminate  the  effects  of  the  diffuser 
in  certain  cases.  This  can  be  accomplished  in  the  present 
analysis  simply  by  setting 


0. 


Electrical  Power  Section 


The  power  conversion  process  in  an  EHD  duct  is  usually 
treated  on  the  basis  either  of  constant  area  or  of  constant 
static  state.  Inasmuch  as  electrical  power  output  per  unit 
mass  is  small,  the  numerical  differences  between  the  results 
computed  by  these  two  methods  is  negligible.  For  definite¬ 
ness  in  this  analysis,  however,  we  assume  constant  area. 

The  electrical  power  that  can  be  obtained  from  an  EHD 
duct  of  constant  area,  negligible  change  of  density  and  optimum 
length  can  be  estimated  from  the  one  dimensional  version  of 
Poisson's  equation  which  governs  the  electrical  field.  The 
solution  is  well  known  and  will  not  be  derived  here;  a  detailed 
derivation  may  be  found  in  Ref.  (3) .  The  essential  result  may 
be  written  in  the  form 


pe  *  7  eEb  *4  v4 

where 

PQ  =  gross  electric  power  output,  watts 

e  =  permittivity  of  medium 
-12 

=  8.854  x  10  farad/m  (for  any  gas) 

E^  =  dielectric  strength  of  medium,  volts/m 
=  area  of  duct  (constant) 

V.  =  velocity 


(10.1 


According  to  the  test  data  analyzed  in  Ref.  4,  the  dielec¬ 


tric  strength  is  well  approximated  by  the  expression 

Eb  =  CB3R3P  4 


(10.2) 


where  CB3  is  a  characteristic  constant  of  the  medium. 


The 


data  show  that  Eq.  (10.2)  applies  to  air  or  steam  up  to  about 
10  atmospheres  pressure  with 

CB3  =  9.49  x  10 3  m2  “K/cmh.  (10.3) 


In  this  study  we  assume  that  Eq.  (10.2)  can  be  extrapolated 
up  to  about  100  atmospheres  pressure.  We  also  assume  that  an 
expression  of  the  same  form  applies  to  other  media  besides  air 
or  steam  but  that  each  medium  has  its  own  characteristic  value 
of  the  breakdown  constant  Cg. 

In  the  two-fluid  system,  the  separate  values  of  the 
breakdown  constants  CBl  and  Cb2  are  usually  known,  but  the 
breakdown  constant  CB3  of  the  resulting  mixture  is  seldom 
known.  In  the  absence  of  adequate  test  data  bearing  on  this 
point,  we  tentatively  assume  that  CB3  can  be  estimated  from 
the  hypothetical  relation 

CB3  =  CxRlSl  +  (1  -  x>r2S2]  dO.4) 


Upon  substituting  Eq.  (10.2)  into  (10.1)and  dividing 
through  by  the  mass  flow  rate,  we  find  the  electrical  work  out¬ 
put  per  unit  mass  of  fluid  in  the  form: 


^.(1^ 


2 

B3 

T 


p4  -  Cd3  (Ts4  -  Ts5  )  (1°*5) 


By  use  of  standard  perfect  gas  relations,  Eq.  (10. 5) can 
be  changed  to  a  more  useful  form.  In  this  connection  set 


<TS4  '  Ts5»  *  iTs  . 


and 


3  = 


£  CB3  Po 


(10.6) 


(10.7) 


where 


PQ  =  ambient  pressure,  N/m 
Tq  =  ambient  temperature,  °K 


Notice  that  3  represents  a  dimensionless  version  of 
the  electrical  breakdown  constant. 

Eq.  10.5  now  yields  the  important  result 


1 


The  corresponding  gross  electrical  work  output  per  unit 
mass  of  fluid  is  then  simply 

w*  =  (10.9) 

Unfortunately,  the  known  value  of  e  and  the  typical 
experimentally  measured  values  of  Cg  are  very  small.  Con¬ 
sequently,  the  dimensionless  breakdown  constant  3  is  ex¬ 
tremely  small  compared  to  unity.  Hence  the  drop  in  stagnation 
temperature  ATg  and  the  gross  electrical  work  per  unit  mass 
W*e  are  also  typically  very  small.  This  is  a  very  basic  and 
serious  limitation  on  the  performance  that  can  be  obtained  from 
an  EHD  generator. 

Treating  the  process  4+5  as  isentropic,  we  may  write 
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(10.12) 


This  fixes  f(M5)  whereupon  M$  may  be  found  in  the 
usual  way  from  Eqs.  (5.16),  (5.17)  and  (5.18).  However,  con¬ 
vergence  should  now  be  very  rapid  since  we  may  take  as  a 
first  approximation 

M5  «  M4  (10.13) 

This  solution  now  suffices  to  fix  all  properties  at 
states  4  and  s4  . 
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11.  Condenser/Separator 

The  primary  fluid  which  leaves  the  condenser/separator 
at  station  6  is  assumed  to  be  compressed  liquid  at  known 
static  pressure  Pg  =  P2  =  P1  and  at  known  static  tempera¬ 
ture  Tg  »  T2  .  The  kinetic  energy  of  the  liquid  at  station 
6  is  assumed  to  be  negligible. 

Temperature  Tg  must  be  equal  to  or  greater  than  the 
ambient  temperature  Tq  in  order  to  satisfy  the  require¬ 
ments  of  a  heat  balance  on  the  condenser/separator.  No 
attempt  is  made  in  this  analysis  to  formulate  this  actual 
heat  balance;  instead  Tg  is  simply  treated  as  a  given  or 
known  quantity.  However,  it  may  be  of  interest  later  to 
investigate  the  effects  of  varying  Tg  on  the  overall  per¬ 
formance  of  the  EHD  system. 

The  secondary  fluid  which  leaves  the  condenser/separator  at 
station  z  consists  primarily  of  a  low  molecular  weight  nonconden¬ 
sing  gas.  It  also  contains  some  secondary  vapor,  but  the  mass 
fraction  of  this  vapor  is  treated  as  negligible  in  this  analysis. 
Thus  the  secondary  fluid  at  station  z  is  treated  as  dry  gas. 

The  thermodynamic  process  in  the  condenser/separator  can  be 
most  conveniently  analyzed  and  modelled  by  employing  the  momentum 
theorem.  We  treat  the  condenser/separator  essentially  as  a  one 
dimensional,  constant  area  device  in  steady  flow.  By  equating 
the  net  applied  pressure  and  drag  forces  to  the  change  of 
momentum  flux  in  the  direction  of  flow,  we  obtain  the  following 
simple  relation,  namely, 
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(p5  -  Pz)A5  -  D  =  (1  -  x)m3Vz  -  m3V5  (11.1) 

The  effective  drag  force  D  in  this  equation  represents  two 
irreversible  momentum  loss  effects.  One  of  these  is  the  loss  of 
momentum  produced  when  the  condensible  vapor  in  the  incoming  flow 
is  condensed  find  brought  to  a  halt  on  the  stationary  surfaces  of 
the  con dense r/separator.  The  other  is  the  additional  friction 
drag  force;  the  latter  is  conveniently  expressed  in  terms  of  an 
empirical  friction  coefficient  cf  .  Thus  we  may  write 

D  *  x  m3V5  +  cf(^p5V52)A5  (11.2) 

Next  we  substitute  Eq.  (11.2)  into  11.1),  bring  all  terms 
to  one  side  of  the  equation,  regroup  terms,  and  divide  through  by 
the  quantity  p5Ag  .  We  also  make  use  of  the  following  auxiliary 
relationships,  namely. 


P5^*5 

p5A5 


P5a5 


y3M5 


and 


m3Vz 


Y2R2Tz 


P5A5  "  Y3M5Mz 


3^3*5 


(11.3) 


(11.4) 


In  this  way  the  following  result  is  finally  obtained  - 


E (x)  »  (1 


Pz 

- 


+  (1  -  x  - 


r")Y  3m5  * 


(1 


-  x>»5 


Y2R2^z 

Y 3*3*5 


0  (11.5) 
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The  quantity  E(x)  serves  here  merely  as  a  convenient  abbre¬ 
viation  for  the  sum  of  terms  shown.  It  may  be  said  to  represent 
an  excess  of  momentum  available  over  that  required  to  satisfy  the 
various  losses  that  occur  across  the  condenser/separator.  In  any 
case,  for  assigned  values  of  all  other  parameters,  we  seek  by 
successive  trials  a  value  of  x  such  that  E(x)  =  0  .  Thus  Eq. 
(11.5)  fixes  the  mass  fraction  x  in  a  manner  which  is  consistent 
with  the  other  specified  input  parameters. 

Of  course,  if  the  various  fixed  input  parameters  are  not 
suitably  chosen,  Eq.  (11.5)  cannot  be  satisfied  for  any  trial 
value  of  x  in  the  range  0  <  x  <  1  ,  and  no  solution  is  physi¬ 
cally  possible  under  these  conditions . 
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12-.  Cycle  Efficiency 


Once  a  solution  has  been  found  such  that  E(x)  =  0  ,  it 
becomes  a  s traightfoward  matter  to  calculate  the  correspond¬ 
ing  overall  cycle  efficiency. 

Because  of  the  pressure  drop  through  the  boiler/super¬ 
heater  ,  the  pressure  =  Pg7  at  boiler  inlet  is  slightly 
higher  than  the  stagnation  pressure  Pg^  of  the  primary 
fluid  at  the  ejector  inlet.  We  assume  that  the  ratio  P  .,/P  7 
is  a  specified  constant.  Hence  with  Pgl  specified,  Pg7 

is  also  known.  Also  the  pressures  ?,  =  P,  are  known. 

b  X 

* 

Consider  the  ideal  pump  work  w^  done  per  unit  mass  of 
primary  fluid.  For  the  present  consider  the  hypothetical  case 
of  a  reversible  pump.  Also  note  that  density  changes  across 
the  pump  are  negligible.  Hence  we  may  write 


* 


w 


P 


(Ps7  -  V 
p6 


(12.1) 


The  ideal  gross  electrical  work  output  per  unit  mass  of 

mixed  fluid  has  earlier  been  established  in  Eq.(10.5).  It  is 

* 

denoted  by  symbol  wg  . 

Let  us  now  denote  the  net  useful  electrical  work  output 

per  unit  mass  of  primary  fluid  by  symbol  wnet  •  It  may  be 
*  * 

related  to  wg  and  w^  in  the  following  way 


w 


net 


Joule/kg 


(12.2) 


where  nx  ,  the  excitation  efficiency,  allows  for  the  small 
electrical  power  expended  to  excite  the  system  and  where  n  , 


the  pump  efficiency,  allows  for  irreversibilities  in  the  pump. 

The  heat  input  in  the  boiler/superheater  per  unit  mass 
of  primary  fluid  may  be  written 

* 

*in.  '  (hsl  -  V  -  <12'3> 

It 

where  hg  is  the  enthalpy  of  the  primary  liquid  at  conditions 
Pg  ,  Tg  as  listed  in  suitable  tables  of  properties. 

The  overall  thermal  efficiency  of  the  cycle  can  now  be 
calculated  from  the  simple  formula 


n 


c 


in 


(12.4) 


Both  nc  and  wnet  are  useful  parameters  which  charac¬ 
terize  the  overall  thermodynamic  performance  of  the  EHD  sys¬ 
tem. 


Eq.  (12.2)  shows  clearly  how  the  ejector  serves  as  a 
kind  of  amplifier  which  increases  the  electrical  work  output 
per  unit  mass  of  primary  fluid.  This  effect  is  shown  by  the 
presence  of  the  parameter  x  in  the  denominator  of  the  first 
term.  Notice  the  beneficial  effect  of  a  low  value  of  x  on 
parameters  wnet  and  nc  .  Unfortunately,  the  value  of  x 
cannot  be  stipulated  independently  in  advance;  it  is  fixed 
by  the  other  specified  input  parameters  as  explained  in  the 
previous  section. 

★ 

Notice  that  owing  to  the  fact  that  wg  is  normally 
very  small,  it  is  possible  in  some  circumstances  for  w  . 
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to  take  on  values  which  are  actually  negative.  When  this 
happens  it  means  that  electrical  power  output  is  insufficient 
to  drive  the  pump. 

The  relationships  developed  in  this  and  the  preceding 
sections  make  it  possible  to  carry  out  systematic  parametric 
studies  of  various  one- fluid  and  two- fluid  EHD  cycles,  with 
realistic  allowances  for  the  various  losses  that  occur.  Such 
studies  can  establish  optimum  design  parameters  and  perfor¬ 
mance  limits  under  various  circumstances.  These  results  in 
turn  can  finally  permit  informed  conclusions  to  be  drawn  con¬ 
cerning  the  ultimate  feasbility  of  this  general  type  of  EHD 
power  generator. 


13.  Calculation  Sequence 

In  this  section  we  summarize  the  various  equations  pre¬ 
viously  derived  in  the  approximate  order  in  which  they  would 
be  used  in  the  calculation  of  system  performance.  Also  listed 
are  the  initial  input  parameters  whose  values  must  be  specified 
in  order  to  start  the  calculation  and  various  further  input 
parameters  whose  values  must  be  specified  in  order  to  proceed 
with  various  successive  stages  of  the  calculation. 

Initial  Input  Data:  Pi  =  P2  =  P6  '  T1  '  Psl  '  Tsi  '  ^hsl  ”  hl)  ' 

R  =  8315  Joule/kg°K 


Y1  "  l1  l"n(Psl7P1)  ) 


(13.1) 


(h 


It 


si 


si 


(13.2) 


(13.3) 


(13 .4) 


Further  Input  Data;  Wv  ,  WG  =  W2  '  T2  =  T6  '  Pv^IV 


CD  = 


wv  vv 

WG  'F2  -  E,vlT6)  1 


(13 .5) 


Verify  that  co<<l 


(13.8) 


Verify  that  M2<1 


Further  Input  Data;  x  (trial  value) 


R3  =  xl^  +  (1  -  x)  R2 
c  =  Y2  R2 

p2  (y2  -  “n 

Cp3  "  x  cpl  +  (1  “  X)  cp2 


(13.9) 


(13.10) 


(13.11) 

(13. 12) 


(1  3.13) 

(13.14) 


(13. 15) 


(13.16) 


(13.17) 


Y3  =  1Cp3  S3> 


s3 


—  [x  Cpl  Tsl  +  (1  *  X)  Cp2  Ts2^ 


'P3 


(Yx  +  1) 


f  (ml) 


=  M, 


1  + 


*1  “  1 


T-  Mi 


2  (y±  -  1) 


(y2  +  D 


f(M2) 


=  M. 


Y  2  “  ^  2 

1  +  -hr-  V 


2(y2  ~  1) 


x  /  s2\  IR1  Tsl  Y2  f(M2> 


i2  (1  -  x)  \psl/  \j*2  Ts2  y1  f(Mi) 


(a1/a2) 


A3  “  (1  +  Aj/A^ 


Further  Input  Data;  PQ  ,  Tc 


S1  -  cpi  in  (v)  ■  Rl  *n  (*T, 


So  -  C_o  in  -  R2  in  (Zsl 


P* 


11  (v/  “  R3  K3  ^  (To  /  "  x  S1  "  (1  "  X)  S2 1 


in  (Psx/PQ) 


P  =  P  e 
sx  o 


(13.13) 

(13. 19) 

(13.20) 

(1  3.21) 

(13.22) 

(1  3.23) 

(1  3.24) 

( 1  3. 25) 


(13.26) 


f(Mx)  «  f 


(1  -  A1/A3)  /P g2\  JR3  Ts3  Y2 


VA3>  f!s2\ 
\Psx) 


lR2  Ts2  y3 


f(M2)  (13.27) 


max  \  y-i  + 


( Y  3  +  1) 

2 


Verify  that  f(M  )  =  f<f  av 

X  IUuA 

For  first  approximation  set  Mx  =  M3 


(13  .28) 


(13  .29) 


(] 3.30) 


Iterate  using  Eqs.  (13-  31),  (13-  32),  (1-3.  33) 


V  *  £  ( ■ 1  +  Mn2) 


(y3  +  1) 
2(y3  -'IT 


F* (Mn)  =  f 


,(y3  -  l) 


(3  -  Y,> 


IY,  ~ 


F(M) 


l(n  +  1)  Mn  F1  (Mn) 


Iteration  converges  to  fix  M 


M.  ./  Y-i  “  1 

g(M.)  =  - \/l  +  (-£- - )  m" 

(1  +  y  Mj2)  V  2  1 


^2  \  Y  2  ”  ^  2 

g(M  )  =  - - - y  \  1  +  (-£, - )  M, 

(l  +  y2  m22  )  \|  2  2 


(13.31) 


(13.32) 


(13.33) 


(13.34) 


(13.35) 
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g(My) 


a  -  i/R3  Ts3  X  l/Rl  Tsl  (1  -  X) 

9  1/  y3  \  g(M1)  y  y x  g(M2) 

(1  -  2y3  g2)  -  \/  l  -  2  (y3  +  l)  g2 


l-Yo(l-2Y,g) 


1  \  ^3  “  ^ 

Eq.  (13.37)  applies  provided  that  g  /  »  —  ) - 2 - 

If  g  =  g®  the  solution  reduces  to 


( Y  ■>  -  1) 


My 


Once  My  is  known  from  Eq.  (13.37)  or  (13.35) 


f(My)  =  My 


L  (XJ) 


( Y  3  +  1) 

Ty3~iT 


f(Mx) 

psy  “  Psx  f (My) 
Further  Input  Data: 


P  -21 
3X  lPsx 


f  (M3)  = 


«  £<M*) 


Verify  that  f(M3)  =  f<ff 


(13.36) 


(13.37) 


(13.38) 


(13.39) 


(1  3.40) 


(13.41) 


(13.42) 


(13.43) 


(13.44) 


For  a  first  approximation  take  =  m 

Iterate  using  Eqs.  (13.31);  (13.32),  (13.33) 
Iteration  converges  to  fix  . 


/^s3\  (Y3  ”  l]  2 

(if]  *  ( 1  +  -V  M3 


T3  *  Ts3/(Ts3/T3> 


P3  *  Ps3'/<Ps3''P3> 


Further  Input  Data;  nD  ,  |  -g 


-  (l  -  nD) 


t3  M3 


Ps4  "  Ps3  e 


£(V  *  £  =teN)fe)f(M3> 


V-S4/V4  /  ' 

Verify  that  f(M4)  *  f<fmax 

For  a  first  approximation  take  M4  *  j  M- 

Iterate  using  Eqs.  (I3.3I),  (13.32),  (13.33) 

Iteration  converges  to  fix  M„  . 

4 


Ts4  =  Ta  3 


(13.45) 


(13.46) 


(13.47) 


(13.48) 


(13.49) 


(13.50) 


(13.51) 


(13.52) 

(13.53) 


(13.54) 


Further 
(CB  «  9. 

CB3  = 


S 

AT 
_ _ s 

Ts4 

Ts5 

Ts4 

Ps5 


l  ♦  M4 


Ts4/(Ts4/T4) 


T  4(Y3  -  X) 
s4 


Ps4/(Ps4/P4) 


-12 


5  x  103  m2  °K/cmb  for  air  or  steam) 


[x  R^  Cg^  +  (1  -  x)  Rg  cb2^ 


e  CB3  Po 


(13.55) 


(13.56) 


(13.57) 

(13.58) 


Input  Data  s  £  =  8.854  x  10  farad/iti,  9  ^B2 


(1  3.59) 


(1  3.60) 


•  8 


Y3  '  "VlsiUj: 


TT- 


Po  i\Ts4 


-1/(Y,  -  1) 


y3  '  2 


1  +  <■  — )  m4 


(1  3.61) 


AT, 


(1  - 

is4 


(1  362) 


P  (-Si) 
Ps4  T  . 


y3/(y3  “  1) 


(13.63) 


f(M5)  -  f  -  / 


/P_„ 

)  f(M4) 

-s4  y*s5  1 

Verify  that  f(M-)  =  f<f 

o  max 

For  a  first  approximation  take 

Iterate  using  Eqs .  d  3* 31)  ,  (1  3.32),  (1  3.33) 

Iteration  converges  to  fix  . 


TS  *  Ts5/(WV 


(v3  -  1) 


Further  Input  Data; 

P  Cf  2 

E(x)  -  (1  -  ~)  +  (1  -  x  -  —  )y3*V 
P5  2  3  5 


(13.64) 

(13.65) 

(13.66) 


(13-67) 

(13.68) 


(13.69) 

(13.70) 


(13 .71) 


Iterate  from  Eq.  (13.14)  until  x  converges. 
Further  Input  Data:  <pgi/pS7>  '  p6  '  ^hsl  ”  h6* 


Psl/(Psl/Ps7) 


Cp3  ATs 


(Ps7  -  V 


(13.72) 

(13.73) 

(13.74) 


net 


qin 


RESULT  (13.75) 


(13.76) 


w 

nc  "  q 


net 


in 


RESULT  (13.77) 
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14.  OPTIMUM  VELOCITY  RATIO 

The  calculation  sequence  suirmiarized  in  the  preceding  sec¬ 
tion  reveals  that  a  coraple  cycle  calculation  requires  us  to 
specify  numerical  values  for  about  thirty  input  parameters. 

One  of  these  is  the  velocity  ratio  (V^/N^)  ~  y  .  The  calcula¬ 
tion  procedure  outlined  above  then  permits  us  to  find  the  cor¬ 
responding  value  of  the  mass  flow  ratio  x  ,  assuming  that 
there  exists  a  value  of  x  which  satisfies  the  governing 
thermodynamic  relations  for  the  specified  values  of  the  other 
input  parameters. 

If  we  now  change  the  value  of  y  but  hold  all  other  in¬ 
put  parameters  fixed,  we  can  repeat  the  above  procedure  and 
find  a  corresponding  new  value  of  x  ,  if  such  a  value  exists. 
It  is  evident,  therefore,  that  x  becomes  some  definite 
function  of  y  as  long  as  we  remain  in  the  domain  where  a 
real  solution  exists.  Of  course  all  other  dependent  variables 
of  the  cycle  are  also  functions  of  y  .  In  particular,  the 
overall  cycle  efficiency  nc  is  some  definite  function  of  y  . 
Moreover,  for  prescribed  values  of  the  other  input  parameters, 
there  will  be  some  definite  value  of  y  ,  let  us  designate  it 
as  the  optimum  value  YQpt  >  which  yields  the  greatest  value 
of  n  that  is  possible  under  the  specified  input  conditions. 
Our  basic  purpose  therefore  is  to  determine  for  any  prescribed 


values  of  the  other  input  parameters,  the  value  of  YQpt  anc*  ^e 

corresponding  value  of  Of  course  the  values  of  x^„. 

c  max  opt 


and  other  dependent  variables  are  also  of  interest. 
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Since  a  working  computer  program  is  now  available  to  es¬ 
tablish  values  of  y  .  ,  )  ,x  .  and  so  on  for  pre¬ 

opt  c  max  opt 

scribed  values  of  the  other  input  parameters,  we  can  proceed  to 
study  the  effects  of  chauiges  in  these  other  parameters  and  to 
search  systematically  for  such  values  of  the  most  important 
input  parameters  as  will  yield  the  best  overall  performance 
of  the  system. 


15.  Conclusions  and  Recommendations 


Based  on  the  work  reported  here  and  in  Ref.  (6) ,  the 
following  conclusions  are  evident: 

1)  Ejector  performance  is  subs tantially  improved 
through  the  use  of  a  high  molecular  weight  pri¬ 
mary  fluid  (i.e.,  above  100)  and  a  low  molecular 
weight  secondary  (i.e.,  below  30).  The  mixing 
penalty  which  ensues  is  offset  by  the  substantial 
decrease  in  ejector  losses. 

2)  Hie  breakdown  strength  limitation  is  severe  even 
at  higher  pressures.  Since  the  breakdown  strength 
of  a  mixture  can  often  be  smaller  than  that  of  a 
single  fluid,  allowance  for  this  must  be  made  when 
considering  fluids  to  be  used  in  the  ejector.  In 
Ref.  (6)  it  is  stated  that  in  mercury /hydrogen  an 
improvement  by  a  factor  of  5  in  the  breakdown 
strength  is  required  for  a  competitive  thermodynamic 
efficiency. 

3)  Consideration  of  the  breakdown  strength  of  mixtures 
together  with  other  factors  seem  to  indicate  that 
the  sequence  of  components  shown  in  Fig.  1.1  might 
not  be  the  best  possible.  In  particular,  serious 
thought  should  be  given  to  an  arrangement  where  the 
electrical  section  precedes  the  ejector.  Here, 
mostly  secondary  fluid  is  present  and  the  breakdown 
strength  will  be  at  its  highest  value. 
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The  termination  of  the  contract  with  DoE  has  left  some 
questions  about  the  two- fluid  cycle  unanswered.  With  the 
analysis  presented  in  this  report  it  should  be  possible  to 
carry  out  further  thermodynamic  studies  with  various  fluid 
combinations.  Also,  more  fundamental  knowledge  on  the 
breakdown  strength  of  dielectric  fluids  will  be  necessary 
before  the  prospects  for  successful  EHD  power  generation 
can  be  fully  evaluated. 


\ 
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16.  Sample  Output 

A  sample  output  is  shown  in  the  following  pages.  The 
case  is  for  mercury /hydrogen  expanding  from  1500  psia  to 
700  psia.  The  value  of  the  breakdown  constant,  CB ,  is  8000 
The  component  loss  parameters  have  numerical  values  corres¬ 
ponding  to  CASE  B  in  Ref.  6. 


FLUID  PROPERTIES 


1 

PS1  =  1500.0  PSTA 

rm=  1^45.0  riFo 

» 

3 

F'l  =  700.0  PSTA 

4 

rii  1770.0  FIFO  R 

5 

(HS1-H1 >=  13.50 

fi  Tl  i/LR 

6 

rw<;i-HA)=  1  70  70 

PfTM/f  Bf 

7 

RH06=  0.750OE  03 

I.R/FT3 

8 

0B1  =  O.ROOOF  04 

M7 

riFR  K/fHR 

9 

NFUJID=  2 

to 

1.10=  O.2016F  01 

11 

GAMMA2=»  1.4000 

1  2 

rp7=  o.rooof  04 

M'"> 

tifr  K/r hr 

SYSTEM  PARAMETERS 

1 

1 

ETAE*  0.85  2 

ETAD=  0.85 

3 

A3/A4=0.?50  4 

CF= 

c 

CD 

c 

5 

ETAX=  0.98  <4  . 

ET AP*  0.85 

7 

PS1/PS7=  1.00 

COMPUTATIONAL  FORMAT 

1 

DELY=  0.10  2 

YMIN=  1.00 

7 

YMAX=  5.00 

WV=  0.2058E  03(F.FF  MOL  UT  VAPOR)  OAMMA 1 =  t.lflIA  Y'*rp=  O.lOM 
ENTER  CONTROL  CODE 


Y 

X 

A1/A3 

OWfDftl  | 

FFF 

REl  Al  fUF 
EFF 

M3 

MS 

5.00 

0.4375 

0.0015 

0 . 0228 

0.7285 

0.0285 

0.0071 

4.90 

0.4531 

0.0017 

0.0213 

0.2175 

0 . 0795 

0 . 0074 

4.80 

0.4375 

0 . 0016 

0.0228 

0.2285 

0 . 0796 

0.0074 

4.70 

0.4375 

0.0016 

0.0228 

0.7285 

0.0303 

0.0076 

4.60 

0.4277 

0.0016 

0.0230 

0.2384 

O . 0307 

0.0077 

4.50 

0.4063 

0.0015 

0.0761 

0,2670 

0.0708 

0.0077 

4.40 

0.4219 

0.0016 

0.0244 

0.7446 

0.0.719 

O . 0080 

4.30 

0.4063 

0.0016 

0.0261 

0.2620 

0.07'>7 

0.0080 

4.20 

0.4063 

0.001 6 

O . 0761 

0.7670 

O . 0770 

O  .  0087 

4.10 

0.3906 

0.0015 

0.0780 

0  .  '1807 

11 . 0  773 

O . 0083 

4.00 

0.3906 

0.0016 

0 . 0780 

O.  >H<>8 

O .0341 

0 . 0085 

3.90 

0.3750 

0.0015 

0.0300 

0.7011 

0.0746 

0.0086 

3.80 

0.3594 

0.001 4 

0.0327 

O .  .V  *  32 

0.0  <50 

0.0088 

3.70 

0.3672 

0.0015 

0.031 1 

0.3119 

0.0167 

0.0091 

3.60 

0.3594 

0.0015 

0 . 0722 

0.3737 

0 . 0  <70 

O . 0090 

3.50 

0.3438 

0.0015 

0.0346 

O  .  <4  7 7 

0.0376 

0 . 0094 

3.40 

0.3477 

0.0015 

0.0340 

0.341  ) 

O . 0708 

O  .  0097 

3.30 

0.3398 

0.0015 

0 . 0357 

0.7537 

0.0 <97 

O  .  0099 

3.20 

0.3281 

0.0015 

0.0372 

0 . 3737 

0 . 0406 

0.0107 

3.10 

0.3281 

0.0015 

0.0372 

0.3737 

0.0419 

0.0105 

3.00 

0.3125 

0.0015 

0.0401 

0.4078 

0.Q4-’8 

0  >0107 

\ 


3.00 

0.3125 

0.0015 

0.0401 

0.4028 

0.0428 

0.0107 

2.90 

0.3125 

0.0015 

0.0401 

0 . 4028 

0.044 < 

0.0111 

2.80 

0.3125 

0.0016 

0.0401 

0.4028 

0.0458 

0.0115 

2.70 

0.2969 

0.0015 

0.0433 

0.4349 

0.0470 

0.0118 

2.60 

0.2969 

0.0016 

0.0433 

0 . 4350 

0.0488 

0.0122 

2.50 

0.2969 

0.0017 

0.0433 

0 . 4350 

0.0508 

0.0127 

2.40 

0.2891 

0.0017 

0.0451 

0.4524 

0.0526 

0.0131 

2.30 

0.2813 

0.0017 

0.0469 

0 . 4708 

0.0546 

0.0136 

2.20 

0.2813 

0.0017 

0.0469 

0.4708 

0.0570 

0.0147 

2.10 

0.2734 

0.0018 

0 . 0488 

0.4902 

0.0594 

0.0148 

2.00 

0.2734 

0.0018 

0.0488 

0 . 4®03 

0.0624 

0.0156 

1.90 

0 . 2656 

0.0019 

0.0509 

0.5109 

0.0653 

0.01 63 

1.80 

0.2656 

0 . 0020 

0.0509 

0.5110 

0 . 0689 

0.01 7? 

1.70 

0.2656 

0.0021 

0 . 0509 

0.5111 

0,0  730 

0.0187 

1.60 

0.2656 

0.0022 

0 . 0509 

0.5112 

0.0775 

0.01<»3 

1 .50 

0.2656 

0.0024 

0 . 0509 

0.5114 

0.0R07 

0.0706 

1.40 

0.2734 

0.0026 

0.0489 

0.4911 

0 . 0890 

0.077’’ 

1.30 

0.2813 

0.0029 

0.0470 

0.4719 

0.0964 

0.0'»40 

1.20 

0.2969 

0.0034 

0.0435 

0.4365 

o.:0$6 

0.0263 

1.10 

0.3125 

0.0040 

0.0403 

0.4046 

0.1165 

O.029O 
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17.  Computer  Program 


In  the  section  the  FORTRAN  program  used  is  listed. 
Results  of  our  calculations  are  given  in  Ref.  6. 


C  T.iU  Piun  £H0  OhO^AT'.U 

IMPLICIT  A  L  *  ■+  (  A-H,U-Z) 

j  )  ,  T  S  l  c- )  ,  A  ,  -  {  6  )  ,  J  1 1  )  ,  T  l  6  )  »  KHL'I  6 )  ,  o( o>  , V l 4, 7)  , 

♦  ,C  d2  .Odd  ,  LF  ,  JELh,  J.iLnlo  .  JELY, 

»  u  ih.  Tl*  a  i  t  TAt-1  » t  T  A.  t  F  X  L,rd2,&_.'lJ,H,''’14+tFA5,FMX, 

♦  >“  3  Y  ,  F  -1 A  X  tr  1  >  F  \  i  G  A'i .  i  A  1  i  jA  ^  'lAil  i  j  *  "1 

+  ■->  1  :j  *  j2  i  *  »  j  23  ,  .2  ♦  » u  2  5  »  j2  u  »  j  j  1 »  j3  2  ,  G3  j  ,  *»3  ♦ ,  j3  3  ,  jJ  'j  ,  j  -l  1 , 0  -1 2  »  G‘1  X, 

+•  <3  a  i  a  ■'  y  ,  ,<K  J  »  a  *\1 ,  *,r-L  J  1  >,?VTo,.»0,*>SlPS7,  jl  ).sl 

♦  :» 3  IX  »  T  J  t  _  i  AlX  ),(kt  1 1\  t:  T  » .i  a  *  .W  »  Y  ,  Y  ■1I'vi,YMiX,T1i»X-12fi)l  »  F  M  Z  r  P  Z  * 

♦  J(  Id  ,12)  ,.,l  2.3) 

C  UA  T  A  F  I  L  i 


J(  1 

,  l)3 

1330. 

Jl  l 

*  2)  3 

13/4.3 

!J(  1 

»  a)  = 

9  00  . 

'Jl  1 

»  4  )  3 

4)1.7 

Jl  l 

,  5)  3 

-♦1.4 

01  1 

,  £)- 

3  76.3 

01  l 

,  /)  3 

♦  /.  17 

01  1 

,3  J  3 

J4oo. 

Jl  1 

,  9)3 

1  . 

Jl  1 

,  10) 

=  13. 

Jl  1 

.  11) 

3  1.3 

Jl  1 

,12) 

39<+90. 

01  2 

,  l)3 

13  JO. 

012 

,  2)  3 

1  Jd9.4 

01  2 

,3 ) 3 

300. 

U(  2 

.4)  = 

ill. 9 

Jl  2 

,5)3 

52-0. 

01  2 

,  J)  3 

73o  .  3 

Jl  2 

,  7)3 

4/.. 35 

Jl  2 

,  3)  3 

9  *t  -)  0  . 

Jl  2 

,  ♦  )  = 

1  . 

Jl  2 

,  1C  ) 

3  13 . 

Ul  2 

,  11 ) 

=  1.0 

01  2 

,  12  ) 

=34 jJ. 

01  i 

,  1)  3 

1 500  . 

01  3 

,2)  3 

i  109 .2 

Jl  j 

,  3)  3 

7  0  0. 

Jl  3 

,4)3 

962.3 

Jl  3 

*  5)  3 

o4.  3 

Jl  i 

» e> ) 3 

7  id .  5 

01  3 

,  7)  = 

>  3.  73 

Jl  3 

,  d)3 

>40  0. 

Jl  3 

,  •»  )  3 

1. 

01  3 

,  10) 

3  l  3  . 

Jl  j 

,11) 

3l.  3 

01  'J 

,12) 

-9<+~i0. 

Jl 

,1)3 

1  iO  j  • 

Jl  '4 

,  <.)* 

i  113.6 

Jl  4 

,3)  3 

4  0  J  • 

Jl  4 

,  4 )  - 

a4p  .  0 

Jl  4 

*  3 ) 3 

7  3.4 

Jl  ♦ 

t  c )  3 

7  7  3.7 

01  4 

,  7 ) 3 

4  4,  73 

JU 

»  d  )  3 

4490. 
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Ul  4,9)  = 
d<  4, 10) 
jt  *♦.  m 
Ul  12  ) 
J  (  n  ,  i  )  = 
Ul  5 , t )  = 
Of  p  *  3  )  = 
Ul  6,4)  = 
Ul  5, 5)  = 
Ul  p*  c )  = 
Jl  5,  7)  = 
Ul  5,fj)=' 
U(  P,5)  = 

Ul 5,  10) 
U<  5, 11): 
d(  p, Li )  ■ 
U  l  O  ,  1  )  = 
Ul  0,2)  =, 
Ul  D,  3  )  a  ' 
U  (  o  *  4  )  = 

ul  5,5)=; 

U(  O,  6)  a  1 

U  <  6 ,  7 )  =  , 
U( »,&)a- 
Jl  6, 9)  =  2 
1J(  p,  10)  = 
U(u,li)= 
Ul 6, 12)  = 
Ul 7, 1 )  =  1 
UC  7*  2) = 1 
U  (  7 ,  3 )  =  : 
J( 7,41=1 
J  (  7 , 3  )  =  l 
J  (  7  ,  o  )  =  1 
U(  7,  7 ) a 7 
’J<  7,  3)=4 
U  (  7 , 9  J  a  2 
U(  7,10)  = 
J(  7, LI  )  = 
Ul  7,12)= 
U(  6 , 1 )  =  1 
U  (  3  f  2 )  =  l 
U(  8.3)=7 
U(  a, 4)  =  i 
Ul j , 5 ) = 1 
U(  0  *  ft  )  =1 
U<  3, 7)  =  7 
•J  (  3  ,  H  )  =  4 
Ul  3, 9) =2 

Ul  J,  10)  = 

U(  3  » 1 1 )  = 
Ul  8 ,12)  = 
U(  *,  l ) a  l 
u(  7,2)=2 
Ul  3 , 3  )  a  0 


a  1  , 

)  )  a  1  3  • 

)  =  1.3 

'  )  =  4‘f  9  0  . 

a  1  50U . 
=1133.5 
=  500. 
=926.7 
=  96.7 
=313.9 
apJ./6 
=9490 • 

=  1. 

)=  18. 

)  =  l  .  3 
)  =  44  90. 
= 1 500 . 

=  1 94p  . 

=  900. 

=  1300. 

=  3. 

=  126.3 
a  7 4J. 
=4000 . 

=  2. 

1=2.016 
)  =  !.*» 

)  a  4000 . 

=  1  500  . 
=19,5. 

=  3  JO. 
=1770. 

=  11. 
=123.7 
=  750. 
=4000. 

=  2. 

1=2.016 
1  =  1.4 
I  =4000. 

1  500. 
I9t5. 
700. 

17  JJ. 
13.5 
130.7 
7  50. 

400 0  . 

2. 

=2.016 
=  1.4 
=40u0. 

1  5  0  0  . 

2  05  0. 
o00  * 
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U(  ^  t  ** )  =  1  7  J  j  • 

L! (  9*51=17.  3 
Ul  ,  «  <j  )  =  i  35  •  3 
0  (  -*»7)  -  75  3. 

U(  *J  ,  -J  )  =  4  0  J  *J  m 
J(  ■*,  ■>)  =2  • 

U(  9* 10=2.015 
U(  9,  U)=1.4 
J(  ■**  121=40  U0« 
J(  I0*ii=15u0. 
J(  Uul^UJ. 
Ul  10 ,3 ) =500 • 

U(  lu  i't)  =  loo  J. 

U  <  1  j  ,  5  j  *  1  '■>  .  5 
U<  I0.fr )*137.5 
Ui  1U»7)=750. 

U (  10*3)*  4000. 
0( 10 ,9) =2. 

J( 10,10) *2.016 
J  (  1 .5  , 1 1 )  =  I  .  4 
U( 10,10=4000. 
y(  l*  i)  =  i. 

V( 1  *  2 ) =1  . 

V(  1 ,  3 1  -  J . 

7< 1*4)=0. 

VI  1  *  5) =1 . 

V(  l,0J*l  . 

VI  1,73*1. 

VI  2, 1) =. 35 
V(  2,2)  =  . 05 
V( 2, i)*.25 
V( 2,4)  =  . 3 
V( 2  *  5 )  =  • ,3 
V(  2  » 6 )  =  •  35 
V( 2, 71*1. 

V( 3, 1)=.7 
V<  3,2)=. 7 
V(  3*  3) *  1  . 

71 1,4) =2. 

V(  j,  6)  =.  95 
VI 3, i)  =  .  7 
V ( 3 , 7 ) *1 • 

V  (  4,  1 )  =  .  9 

V  (  4 , 2 )  *  1  . 

V( 4,3) =1. 

V( 4,4) *0. 
V(4*5)=l. 

V  (  4  *  o  j  =  .  5 
V14,  7)  =  l. 
w( 1,1)=. 5 

a(  1 , 2 )  =  1  . 

W( 1, 0)=2. 

W<  2,  1) =.  1 
rti  2,  c.  )  =  l  . 

0( 2 , 3 ) =5 . 


c 


.DATa  INPUT 
50  0  t\AlTc(5,5l3l 

510  riJKiAl(,J,»GA,*Eifct'  C  uNT  k  0L  C  uD;  *  ) 

K  c  i  J  (  5  t  p  J  i  i  J|1  Vil* 

52  0  FORMAT { 3121 

IFlIU.t'J.JI  ST  JP 
IF ( IU.c  1.-2  1  GO  TJ  133 
IF  (  IU.iJ.-3  )  GO  T'3  .300 
IF  I  IU.  =  jJ  TO  232 

525  wklTuU, 530) 

5.30  FO  -MAT(  •  0*  ,jX,*E‘iT£R  FLUID  PROPERTY  CHANGES,  IF 
■<&AO(5*54ji  J'JtiU 
540  F  0  k  '  1 3 1 ( I2»F9.4) 

IF  ( J'j.cO.JJ  GO  TO  550 
tFUU.r3.-l»  GO  TJ  500 
IF  (  JJ.ti  j.-2  )  GJ  TJ  133 
IF ( JU.cv.-3 )  Go  Tj  303 
IFUU.C.1.-4)  GO  TO  232 
U(  IU  ,JJ) =  2U 
PU  rc  525 
550  nK I Tc ( Ot  5aO 1 

5t>0  rja-1AT{*0*,3X*»£rjTs:R  SYSTEM  PARAMETER  CHANGES* 


RUAO(5t540)  JVtZV 

TUI  HJ  :  •  .11  .in  T 


IF  (  UV 

01 

GO 

T  3 

5  70 

IF  ( JV 

•  2  j* 

-  L  ) 

Go 

TO 

500 

I F  (  J  V 

-d 

GO 

TJ 

130 

IF  <  JV 

•  il  -J  • 

-3  ) 

GO 

T  J 

303 

IF  (JV 

•  l.;. 

-41 

GO 

TO 

232 

V  (  IV, 

J  J 1 *Zv 

GM  TO  553 
57  J  /.RITEU,  5oJ) 

530  FJkOAT  i  •  3*  ♦  3-<  *  •  E4T2K  FORMAT  CHANGES,  IF  ANY.1  ) 
*  b  AO  (  5  ,  5-»U  1  J  a  •  Z 
IF  (  Jh.c.j  .3  )  GO  TO  130 
Ir(  J.-..c  J.-U  GO  TJ  503 
IF  (Jw.cG. -21  Gu  T<J  130 
IFlJw.cw.— 31  Go  TG  300 
IF  (  Oh.  E'j.— 4  1  GO  TO  232 
w(  Iw  » J 1  -La 
GO  TC  570 
130  PS(  1  1  =  0 (  IJ»11 
T  S ( l >  =  J(  IU  f 2 ) 

PI  1 ) -0 (  I J »  3  1 
I  (  1 1 »U(  I  0 » 4 1 
GcL.HaUl  U,5) 

J£  L  H  i  o  =  J  7 1 0 » o  1 
i50n(  6)=’j{  I  J,  7) 

G j 1=0 (  10,  J) 

GA  •  1AI  =  l .- ALGG (T Sill /Till  1  /  ALJG(?S(  l)/,(l)l 
:M A  1*  1  •  /oAMMA  1 
G 1 1  -  (  Ga.'p  Ma  l - 1 .  1  /  2  . 

G 12s oil f  1, 

Gl  J=ijA  ;i  M  A  1  /  (  G  A  MM  A  L-  l  .  1 
u  1 1*  G 1 3  —  1  • 

G 15=01 2/ ( GAMMA 1-1. 1 


ANY.  •  1 


IF  ANY • ' 1 
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CJ  l=L)tLi/i  i  5  (  i  )  -  1  (  i  )  ) 
iU  =  CPl/  i  1  3 
<•  V  =  l  .  *-  3  5  «  7/  .<  I 
*FLOiu=ii-l  X(  G(  IJ  ,0)  ) 

Ir(  >!H  Jl  J.-o.ll  v,J  TO  4G0 
<•  j  =  ,J  l  1  J  ,  10  ) 
v,.VMAc=u(  1'J.lll 
:52=J( i ,,12) 

GO  TG  '»10 
40 o  Cj2~C31 

G.\y  M6a=OA-MAI 
-  u  V 

410  E  T  Ah  =  V(  I  V,  l  ) 

J A  J—  V  t  I  V, 2  ) 

A  i  A~»=V  (  I  V  t  3  ) 

C  p  s  V  (  I  V  »  4  I 
Er*X=V( IV, 51 
E  T  P  =  V  l  17,6) 

?3LPS7»  V(  w.  n 
Oc  LY  =  .%(  I I  i 
r"U\  =  h{  I  -,,2) 

Y  -1 A  X  =  v>  (  1**3) 

PJ=l4.o9o 
r 0=520. 

1 T  c:  (G,U0)  Pod)  .  IS  t  I)  ,°l  1)  ,T(  l)  ,  GELH,  JELHU,  KrtOl  6)  .Colt 
♦  Jr  LU  1 0  ,  at  J  ,  o  A  ••  1A  2  *  CO  2 
l-,G  FJa  •  J*  ,  3<,  *FL  j;  J  P  31J  P  c.  K  T  I  E  3  •  /  / 

I  PS  1 =*  » F7 , 1 , •  PSIA'iLOX,'  2  T  S 1  =  '  ,  F  7.  I,  •  DEG  R  1  / 


^  P  i  =  •  » r  7 .  1 1 '  PSlP'ti^A,’  ♦ 

+  4X,'  5  ( .IS  L-rti  )  =  ' ,  F  7.  2  ,  *  '->TU/LJ'»S 

otu/lo'/ 

+  4X  , '  7  H-huc-*  ,E  12.4,*  LB/rfONSX,' 

♦3*// 


PS1  A'  ,  hA 
r.2,'  6  TO /L; 


(  iSl-Ho)=* ,f 7.2, 


.M2  L)EG  K/ C>i 


MFL'JIO*'  ,  12,2 IX, *10  av=«  » E 1 2 .4/ 

GAM  3A  2= '  *p [ . 4 » Iox »  '12  w52  =  ,,El2.4,' 


M2  DEG  K/CMS*  J 


rt4ITt(6,  1 5  u  )  E  T  A  r  »  E  T  «u,-»3A4,CF,FrAX,ETAP  ,  P  S  L  P  S  7 
150  Fja.MATI  '  j*  ,3X,'SYsrF  i  P  3  KA  Mr  r  E  .<3 '  /  / 

♦  4X  ,  1  1  ET  Ac=  '  ,F  ->.2,  iX,  •  2  u  T  Ai)  = '  ,  F  5.  2 , 3X 
H  3  A  j  /  A  ♦=  *  »F5. 3,2 X,'  ■+  CF=  '  ,F 3. 2/-»X,  •  5  £rAX=',F5.. 

+•  6  Fr\P= • ,F5.2 , JX, »  7  PSl/P j7*' ,F5.2I 
>.<<  IT2(  j,  153 )  OEl>  ,  f  Y  l  N ,  Y  '1 A  X 
155  F  JKMAT  {  •  O'  ,  3X  , 'C-.lPuTATIu  ,AL  FOR-nAT*// 


♦•AX,'  i  JcL  i- '  ,  F  ”j  •  2  ,  jX  , 
150  Y=Y41N 

p ( 2 ) =p  i  n 

P(61-P12I 
P  2  =  l«c35S7/li)G 
T (  2)  =11  LI 
ri6i*riu 

G2  i=(GA  14A2-i. 1/2. 

22?.=  G2  1  *1 • 

G2j=  GA  Ma>/UAMMA2-U» 
G2 t=  o 2  j  ""  1  • 

32  5=  G22/CUAMMA2-1. I 


Y  M  I M  = '  ,F5.2,  jX,' 


Y 1AX  =  * , F  5 . 2  I 


o26=G25-l . 

A  1lli  =  l T  Si  1 ) /T{ l J-L. 1/Gil 
x  k  i  )  =  i  ;*t  i:<  i(  i) ) 

*2C'-.=x  HU  ■»■«-■!  T  Al./ GA-1  '*.\2*  HL/iiZ*  T  (  1  )/T  (21  I 

r-  u=x;-(  U’H  u  >**<  -til  3 » 

j.'U=*M  U/  i  1J**2)*33xT(  l.«-31L*XM(  1)  **Zi 

<IK  I  T“  {  >iV  i  j  A  X  4m  1  t  Y  2C  K 

230  Fjh  1«T{  •  0*  ,4X,'.tV='  ,c'12.4,  MtFF  MOL  WT  VAPOR!  JArtMAl*' » F  7 .4  «  ’ 

♦  C8  - 1  »F  7  • 

jU  To  500 
MOM  (i. \i  rU'-f  FUNCTIONS 
800  Y  =  Y MAX 
til  -j  ak  l  T t;  (  .» *  32  J)  V 
62 0  F0k.-UT(  •  3*  ,<*X,  *Y  =  'f  F7.3) 
x=  o. 

83 J  JO  830  LOOP=lf21 
XSfir  1 s  c  X  0  fc  S  S 1  X ) 

WaI  TE  (o  .  i40)  AtX  SMTM 

84u  FCF  iAT(  3A, 1  <=  •  ,F7.3,  '  XSM T M* • , 2 1 2. 4) 
x=>.  35 

850  CG.\riNj6 

lr  C  i  .LE.  YMlfsl)  GO  TO  500 

Y= Y-DtLY 
GO  T(j  -J  1 J 

COvJ  IlIsJjUJ  1  AftO  PErlPOKMANCfc 
2U  *RITE( 8,235) 

235  FJP.-1AT  l  •  O'  ,3<,  •¥',  iX, «  X'  ,5X,  'A1/A3'  ,3X,  •  OVERALL'  r 

♦  2X, '  PElATI  v  =  *  ,3X  » •  13*  ,<JX,  »M5V 
«-37X,  *ErF  * , 7X, ' EFF • ) 

Y=Y MAX 

2iG  >=X2£FU(Y) 

30  J  PS7=P5< 1 J /FjI PS7 
tt:*CP3- JcLl  3 

*F=l*t4./778.26*( PS7-P1  6) ) /RH01  6) 
r;EX=«t/X 

IF  ( XLX . Lc. JtCH )  GO  TO  340 
w£  X  =  Oc  LH 
Apire  (  6 1  32  J  J 

320  F 0 k  -1 A T  1  •  j'  ,*X,»  AVAILABLE  ENTHALPY  D^OP  ATTAINED'  ) 
i**’j  a  l£  T  ®£  TA  X*'.vE  X—iaP /£  TAP 
2IN=0hLHlo-wlP/ET  AP 
cT  AC ■=*..'!£  T /  jl  1 
a.-I  A  A  =  t  T  A  .<  *  0  c  L  H—  *  P  /  E  T  AP 
•1 7 AK  =  ir.  It  T/ a  TAX 

akITF  Y,A,U43#ETAC.ETXk,XH3)tXrf(5) 

35 J  FJ8MAT1 *0' ,F0.2, orl0.4> 

Y= Y-OFLY 

IF( Y.Gt. YTIN)  TO  238 

GO  TC.  500 

END 

FJXCTICN  XMACh 

FJNOTI.iN  a-:ACH(F  1N.GI 
ILLICIT  8cAL*41  ^r.l,  u-Zl 
MFL AG=0 
XOUNT* J 


ul®<G-i.)/2. 
o2— 'j1  +  1. 
j5  =  G2  /  (  J-l  .  J 
—  3—  l  . 

XMM*0. 

r  UX=Gi**(-G5 ) 
r-1urvXsF  -IN/ F  MAX 
IHF^r  IX.Ll.i.l  JO  TO  1300 
X  '1  ACri®  1  . 

nMIT-I'j.  13U)  F-1NFMX 

1010  F  1PVAT  l  •  )•  ,3X,  •  F.-.3/F  'UX»'  »F7 .5  ,•  CHOKING  LIMIT  SXCEEJEJ') 
-ETU^K 

1 30  J  F»  F’-VX*  i  L  •  ♦  l»  L  * ^  *•  X  MN 

if  (  vfla  ;.to. u  .  ..Ki t:< ^..iojo)  xmn.f 

luiJ  FORMAT (  '  J'  ,d<,M  =  ',F  1.5,  '  F=',ci2.4) 

JtJ-F^AaS  <  FJ 

IFIJEIF.LT.. jF— 0 o )  GO  TO  1320 

F3=F,-ii\Vo2-X  i  !~(  1.  **36-  l. 

X  4\=  XM.-F/FP 
<JUM=.\GU  4l>i 

[ F ( KlUoT . L>-  «  2  J )  33  TO  1300 

VUTFlo,  l  WJ)  XMW.F 

10*3  F  JrtMATi  '  O'  ,4a,  •  J  =  SF  7.  5,  •  F=',tl2.4,»  AFTEa  20  I  TER 47 1  CMS  *  ) 

1  020  xMACri=X  in 


ritT'Jnm 

FnQ 

C  Fu.vCTION  X2ERJ(Y) 

Fjr.'CTlJn  XZEXOm 
IMPLICIT  rsCAL*4  (A-H.O-Z1 
I T  EA=0 
AR  =  1  . 

E=  cXCbj  S {  l  .  ) 

A3Sc*A33(E) 

IFIAOSE.GE. .5E-04)  30  TO  TOO 

XZ  c*C* 1 • 

RETURN 

733  IFIE.GT.O.J  GO  TG  723 
«XITb(5»710)  L 

713  FJkm.UC  »  J‘ ,4X,«X*l.O  E* '  »  E12 . 4.  * 
XZEROL. 

*E  T J*N 
723  XL  =0 . 

u--XcciS(3.  ) 

AtJS==AJi.(c  ) 

IF ( AUat.Jc. .5E-04)  GO  TO  730 
X  ZERO* -3. 
aE  TUimN 

/‘JO  IFIE.LT.J.)  GO  TO  750 
/.il  I T  c  ( o  »  74  j  )  E 

740  FUk.MATI  4  34 ,4A,  *A®0.0  E  = 4 , El 2 . -*» 4 
aZEiO's  J. 

RETURN 

753  x*(XL*\1)/2. 

E®  EXLEi S ( X ) 

A3SE®A jS ( E ) 


MOMENTUM  JErlCIT*) 


MO.-IEN TUM  I N  VE RS  I JN » ) 
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r- 


7-)  j  r; 


1  r  nP  =  i  ts  xn 

IF  l  I  T  L.. .  1_L .  13  ) 
«./ 1  TF  ( o  *  l  i  0 )  A ,  E 


/.AT  (  *  J*  , 4a , ' X ■ 


JO  TO  7RO 

F  7.5,  • 


c= '  *612.4, 


AFTER  13  ITERATIONS') 


N 


X2EXOA 
ac TURN 

7du  i r ( A /S 2 .GE . • 5£-0 4 )  JO  TO  760 
XZENC-A 
;\E  TURN 

76J  IFlt.LT.G.)  GU  TO  770 

AR  =  X 

•JG  TO  750 
770  XC-X 

GO  TO  750 
SNO 

C  rJ.NCTICN  EXCESS!  X  I 

FJNCT1  J  J  t  aCESS I  X  ) 

ILLICIT  ^c  At.  *  4  (  A  -H  ,  0  - 1  ) 

COM/PN  Pit  &  )  .TSI  o>  ,X  Mo)  ,?(o  ) ,  TI  a)  ,?.H0<  o  )»S(  o )  ,Vl  4, 7)  , 
f  4LA3,A3-w,CPi,2i2,CP3,2Pl  ,232  *  w  3  3  ,  C  F  ,JcLh,0wLHI6,0ELY, 

+  uiLTS,jdLTS-t*~,cTAE,:';TAO,cTA<*'iTA'>,ETAC,r'-;l,F/0»r33»F/4,F'-15,F/Xt 
f  F  /  Y  ,  F  i*>  A  A  *  F  .1 3  F  /  X »  F  /  i  r  ".A  ,  G  A‘i  .1 A 1 ,  j  A  Vi  \  0  ,  i>  3  .‘1A  3  ,  1,->12,j13,J14,o15, 

f  G  16  *  C 2  L  ,  'il e. , 03  ,  o  2 1 »  32 6 » G2  6  *  03  i. » 032  *  Gj>  ./»  G3*t »  J3  5 » uJ  o  » G.'*i  *  6'32  » GMX , 
F  A  /  X  ,  A  /  Y  ,XXu*aI*>1,  i  F  L  J  I  0  t  F  V  T6,P.)»P3l.P57,  }  I N  •  r.  I  ,  i\  2  ,  R  3  , 
Fo3iX»li)»*c***cX»WP»“ti'ifcf  «  aG  »  wv  »Y  *  7  / 1  .v  »  ^  M 1 X  ,  T  2  *  XM  2 »  0 1  iF^ZiFI, 

+  J(  10 ,12)  ,A< 2,3) 

\.i  I  2  )•=  I  A  U  L  )  /  Y )  *  S.jR  T  t  (  J  A  VIA  1/  GA.T1A2  )*IRl/R2)MTll)/T(2))) 

IF  (XM(  2)  .27.  U  )  JJ  TO  234 
a  -I  C  2  )  =  1 . 

1  ( 2  Nl'U  )/Y**2*JAM •( A 1  /  GAM -IA2 * f\l/R  2 
23  i  TS2T2=L.+J2L-XM( 21**2 
ra<2)=  TIZJ-'Ti^rj 
Pi(2)=  P(2)*TS2T2**G23 

G/2-A.V  I  2  )  /  (  I .  +-GA.  !MA2  *X/I  2)  **2)  ♦SORT l  1.  «-G21*XM(  2 )  **2 ) 

240  K3*A*r4  +  (l.-XI*i'2 
CP  2= J2 J*R2 

CP3=a*CP1M  l.-X) *C?2 
G  A/M  A3 =CP3/ (CP3-R3) 

G3l=(GA.MA3-U)/2. 

J3  2-G3  LfI . 

G33  =  GA.HA3/(  JAM 4*3-1  .  ) 

G  3  4- Gd j- L . 

G35-G32/ (GA 3 tA3-i.) 

Jio=J3S- L* 

F  -1AX  =  'j j2*  s  (  -G35 ) 

TS(3  )  =  i  <*CPl*T$(  1  )  +  (  L.-X)  *CP2*TS(  2  )  >/C.°  3 
F.4c=x:*:(  2)  M  1 .  ♦  J2  1  *aM  I  2) **2 ) ** I  -J2  5  ) 

3l  =  PS(2  )  /  Pi  U)-*F  v2/F/l>f‘S:a  Tt  R1/R2YTSI  1)  / TS I  2 ) *GA  4  4A2/GA 4/A i ) 
AlAJ»X*Jl/I  1.-Xfx*01  ) 

SI  1 )  »CP  l  * ,  LOGITS  »  l)  /  rj)-Rl*AL-J(PS(  D/^O) 

SI  ?.)  =  C?2  *AcOG  I  T  S  (2)/T0)-.42*ALjG(P3I2)/PO) 

PSX*(CPi*ALJa<T.>l  jj/T0)-X*Sl  l> l  .-X ) *S< 2) )/ X3 
PSXaPO-h  MI  PSX) 

02=Pit2)  /PiX*F/2*SJi< T(R.3//2*TS(3)/TS(2)*GA.-i.|42/GA/^A3) 

F4X*D2/I I.-XFX*Jl) 
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5  ,  *  FMX/FMAX*'  ,  £12.41 
QRT(  <l*TS(  1)  /GAM  AAD+l  1.- 


ir (Fvx.ie.r 4axi  3a  to  2vi 

X'1X  =  L. 

F  -1 X H  =  F X a / F  *AA 
F  4X  =  FMAX 

1  T  t  (r.  .243)  Y»  )  iF  i^XFMX 
24j  FU.<-:.»T  {  *  J*  ,4A, '  V*'  ,F  7.5,  *  X  =  ',F7. 

GJ  TC  2  4- 

241  X  -1 X*  XM  -  G  H  (  FM X  ,  G A*-  4  A  3  ) 

24v  G.^Y»Sw.;r(.n*ri(i)/aA.vYA3»/  (x/gu*s 

+  G  jkT  (K2^T..(  2  I  /(,«;•.  1A2J  ) 

g--iax*u/sjk t  u.*u+:'  -ia3*i.  n 

iP  (C.'1Y.Lc.o4AX)  ra  2 45 

J«1Y  jVXsG.'iY/G  -MX 
n  R  I  T  c  <■.»,' 42)  Y  ,  X  ,  GFiY'GMX 

2-2  FORMAT!*-)'  ,4A,»Y  =  *,F7.5,*  X=*,F7.5»*  GMY/GMAX**  ,  512.4/ 

♦  5X  t  *  ."G  1e..)T0f'  LIMIT  EXCEEDED. •) 

X'<Y*1. 

f-  M  Y=Ff*.A  X 
GJ  TC  24  3 

245  X  !  Y  -  {  (  1  .-2. *GA«MA3*GMY**2  J  -SORT!  1 

*  (  1  ,-uA  IX  \3<M  1  .-2.-J  m  HA3*"jMY**2  I  ) 
amvsS  ;,kT  (  <  -5  0 

FlY-x  '1/(1  .  -.,3l*AMY**2i**G35 
IF  (F Tr.LE.r.-UXJ  oJ  TO  243 
X  ’•*  Y  —  i  « 

rJlrnspNI/FlUX 
F.  1  r= Fii  X 

W*ITc  i  6.246)  Y.X.FM  YF  3X 


-2.*l  GAMMAa+U  )  * GMY**2)  )/ 


2-6 
24  j 


F  /  .  5 , 


f  1AT{  *  J* ,4X, 1 Y-'  t 
p>SVsPS  A^FlA/r.-tY 

{  i)*r>G  <*(  (  PSY/PSX)  **(  1./  ETAE  )) 

F  IJsASn/^SI  3»*FMX 
IF  IFV3.LE. F1UI  GJ  TO  249 
X  U  3  I  *  1 . 

,yX  I  T  L I  j  t  3  JO  J  X  4  (  11  |X  *(2)  t  <  1X.XMY.X 
♦  F.UX  •  .i .4 2  » GMY  i  u’tuX t  PS  (1)  » P  S  1 2  ) » 

*TS(2),T3(3) 


X  =  *  tF  7.  5  »  *  FMY/F.-1X**  ,S12.4) 


X(  3)  ,F -11tFM2  ,r  -1X,F^Y,F^3t 
PSX»FSYi>*S(3)tTSUIi 


4vKJ  F'.IX'UT  (  »  J*  ,i>X  ,Mi,  )2.MX,MY,M3-*/5512.4/ 

**D*  f  *F  11  ,F*,2»F4X  ,FMY  .  F  ,-U  ,  F  1  AX  =  ’  /  5  E 1 2  •  4 »  E  1 1 .  -/ 

♦  *  o  *  t  *  G.-i  1 1  »'•)  2  *  vj  T  Y  t  J  1 A  X3  *  '/4E  12 . 4/ 

•f*  3'  r  *  °3l  ,Fi2,ySX  .PjY  ,PSi=  *  /5 6  12.4/ 

♦  *J'»*TjL,rj»2.TS3  =  */3el2.«*) 

F  I  J  F  '**r.4)/F'1AX 

F  ‘4  3  a  r  Pim  X 

.RITE  (j.'il)  Y  t  X  » F  M  3  F  M  X 
25a  FCinMaT  (  •  J*  t  4A  t  *  Y*  * » F  7. 5,  *  X  =  »  ,F7. 

GO  TO  2*  7 
249  X,t(  3I*a  na.-i(F.-}3#GA:-u4A3J 
2-7  rSJTJ*L.*GU*X  U  3)  **2 
T ( 3 1  —  T i ( 31/ TS3I3 
PSJP  i=TS  iTj*?333 
P(3)*F5»3)/PS  J  P  3  • 

Pit  V)*>>  »{  3)  *2API -(  L.-tTAO)  *GAXMA3/2.*XM(  3)**2) 
F'iVsFSl  3) /PS  I  t)*AJ44*rMJ 
ir  ( F M4  .  Lc  «F  TAX I  G  J  Tu  252 


5,'  F  V3/FMAa=  *  *  El 2.  4) 


XJ/GM2* 
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2=53 


252 


25-t 


X  '3  { 4  )  =  1 . 

F.44F*X*|M4/F;3AX 
F  .34*  F  ("AX 

**17P<„,253)  Y  f  X  »  F  14FHX  „  . 

F'J.<‘‘!4  H  '  J*  »  ♦.<  ,  *  Y*'  ,  F  7 . 5,  *  X*'  fr7.i»  1  FM4/F  3A  X  =  '  ,  £12.4) 
j»J  Ti J  2  54 

X  3  (  4  )  *  X'3  A-  H  (  F  H4  ,  0  A.3M  A3  ) 

73 (4 ) =7 313) 

7  3  47  4*  L  «  +0  3  l  vXi*l(  4  I  **2 
7 { 4) *T 3( 4 ) / 73474 
Pi-«P  ~i—  To  4T 
P(4)=P:>(4)/  !*  34P4 


Cjj=X^iUU  l.-X)*Cd2 

TA*0.  JGJ9o75*<  pJ/14.69o)  *<CB3/  7490. *r>2  O./TjI  **2 
0EL7S4*.>t7**G  U/OA  )Y..3*?S<  4)  /p  0*(  TO/ TSl  4 ) )  **  2/  11.  +  G3 


1*X:M4)**2  )** 


255 

253 
2b  D 


+•03  4 

0ELTS»ril4)*JELTS4 
73 (5 )  =  73 14)—J=L7 5 
P3(3)*P3(  »)*<  Ti(  3  )/l  3<4)  )**G33 
F35=  fit  j  )/7i(  4)  «=PS(  4)/PSI  5)*FM4 
IF l  FM5. 1_ E . rA\K )  00  TO  250 

X  '3  (  5  )  =  1 . 

F.3  3FivX*F:35/F,3AX 
F.35  =  F.v.aX 

W, UTE  (0,255)  Y ,  X , F M5F'3X 

FU.s'-iAT  l 1  C»  ,4<f  »  Y=*  , F7.5,  •  X='  ,F7.  5,  ' 

00  TC  260 

X. U5)*X*»ACH  (  F/i5  »0AMKA3) 

TS5T5*  l  .  ♦  jjlfA.q  5  )  **2 

T  (  3 )  =  i  0  (  5 )  /  T  05  75 

P  5  5P  5=  T  0  5  7  b*  +G  3 3 

P( ?)*PS15)/P0pF5 

F3Z=l l.-X) *AiA4*FM2/< i.-X+X*Dl ) 

a.32=Xi’1mCH  (  F  32  1  Ga.v'’3a2  ) 

7S277=1.+G21*XH2**2 


7Z=7i(2)/7S272 


FM5/F-UX** 


PZ  =  ?3(2)/U3Z7Z**023)  „  B 

£XC£SS*L.-P2/P(pi  +  (  l  .-X-CF/2. )  «»A  3?*43*XM  ( 5  )  **2-(  1 
+  X.K  5)*X.3Z*3  jn7(OAM3A2/0AM  3 A  :>** 2/ * 3 *  7 Z/T l  5)  ) 


RE7UKN 


£.30 


£12.4) 


-X  ) <=GAMMA3* 
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18.  Nomenclature 


A 

— 

area 

a 

- 

speed  of  sound 

CB 

- 

characteristic  breakdown  i 

S 

- 

specific  heat  at  constant 

cf 

- 

friction  loss  coefficient 

Eb 

- 

breakdown  electric  field 

E  (x) 

- 

Eqn.  11.5 

F(Mn) 

- 

Eqn.  5.16 

f  (M) 

- 

Eqn.  5.3 

^max 

- 

Eqn.  5.15 

g(M) 

- 

Eqn .  7.5 

g 

- 

Eqn.  7.7 

gmax 

- 

Eqn.  7.9 

g- 

- 

Eqn.  7.10 

h 

- 

enthalpy 

I 

- 

irreversibility 

M 

- 

Mach  number 

itl 

- 

mass  flow  rate 

Pe 

- 

gross  electrical  power 

P 

- 

pressure 

^in 

- 

heat  input 

R 

- 

gas  constant 

R 

- 

universal  gas  constant 

s 

- 

entropy 

sm 

- 

entropy  of  mixing 

73 


T 

V 

W 


w 


w 


w 


net 


y 

2 


temperature 

velocity 

molecular  weight 

gross  work  out  per  unit  mass 

ideal  pump  work  per  unit  mass 

useful  electrical  output  per  unit  mass 

*1/A3 

vl/v3 

quality  of  vapor 


8  -  Eqn.  10.7,  dimensionless  breakdown  constant 

r  -  ratio  of  specific  heats 

A(  )  -  change  of 

e  -  permittivity  of  the  medium 

nc  -  overall  cycle  efficiency 

nE  -  Eqn.  8.13,  ejector  effectiveness 

nQ  -  Eqn.  9.6,  diffuser  effectiveness 

n  -  pump  efficiency 

1? 

nx  -  excitation  efficiency 

p  -  gas  density 

!|»  -  steady  flow  availability 

u  -  mass  ratio  of  condensible  primary  to  dry  secondary 
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m 


Tl 


T 


Subscripts 

0  -  ambient  conditions 

1  -  primary,  entrance  to  ejector  (exit  from  boiler/superheater) 

2  -  secondary,  entrance  to  ejector  (exit  from  separator) 

3  -  mixture,  exit  from  ejector  entrance  to  diffuser 

4  -  mixture,  exit  from  diffuser  entrance  to  generator 

5  -  mixture,  exit  from  generator  entrance  to  condenser 

6  -  primary,  exit  from  condenser  entrance  to  pump 

7  -  primary,  exit  from  pump  entrance  to  boiler 

f  -  saturated  liquid 

G  -  gas  or  secondary 
g  -  saturated  vapor 
n  -  next  value 
s  -  stagnation 
v  -  vapor 

x  -  exit  station  of  hypothetical  minimum  loss  device 
y  -  exit  station  of  ideal  ejector 
z  -  exit  station  of  condenser 
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